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Abstract 

We reformulate the analysis of nuclear parity-violation (PV) within the frame- 
work of effective field theory (EFT). To 0{Q), the PV nucleon-nucleon {NN) inter- 
action depends on five a priori unknown constants that parameterize the leading- 
order, short-range four-nucleon operators. When pions are included as explicit 
degrees of freedom, the potential contains additional medium- and long-range com- 
ponents parameterized by PV ttNN coupling. We derive the form of the corre- 
sponding one- and two-pion-exchange potentials. We apply these considerations to 
a set of existing and prospective PV few-body measurements that may be used to 
determine the five independent low-energy constants relevant to the pionless EFT 
and the additional constants associated with dynamical pions. We also discuss the 
relationship between the conventional meson-exchange framework and the EFT for- 
mulation, and argue that the latter provides a more general and systematic basis 
for analyzing nuclear PV. 



1 Introduction 



The cornerstone of traditional nuclear physics is the study of nuclear forces and, over the 
years, phenomenological forms of the nuclear potential have become increasingly sophis- 
ticated. In the nucleon-nucleon (NN) system, where data abound, the present state of 
the art is indicated, for example, by phenomenological potentials such as AV18 that are 
able to fit phase shifts in the energy region from threshold to 350 MeV in terms of ~ 40 
parameters. Progress has been made in the description of few-nucleon systems [1], but 
such a purely phenomenological approach is less efficient in dealing with the components 
of the nuclear interaction that are not constrained by NN data. At the same time, in 
recent years a new technique — effective field theory (EFT) — has been used in order 
to attack this problem by exploiting the symmetries of QCD [2]. In this approach the 
nuclear interaction is separated into long- and short- distance components. In its origi- 
nal formulation [3], designed for processes with typical momenta comparable to the pion 
mass, Q ~ m n , the long-distance component is described fully quantum mechanically in 
terms of pion exchange, while the short- distance piece is described in terms of a num- 
ber of phenomenologically-determined contact couplings. The resulting potential [4, 5] 
is approaching [6, 7] the degree of accuracy of purely-phenomeno logical potentials. Even 
higher precision can be achieved at lower momenta, where all interactions can be taken 
as short ranged, as has been demonstrated not only in the NN system [8, 9], but also 
in the three-nucleon system [10, 11]. Precise (~ 1%) values have been generated also for 
low-energy, astro-physically-important cross sections of reactions such asn + p— * d + 7 
[12]. Besides providing reliable values for such quantities, the use of EFT techniques 
allows for a realistic estimation of the size of possible corrections. 

Over the past nearly half century there has also developed a series of measurements 
attempting to illuminate the parity- violating (PV) nuclear interaction. Indeed the first 
experimental paper was that of Tanner in 1957 [13], shortly after the experimental confir- 
mation of parity violation in nuclear beta decay by Wu et al. [14] . Following the seminal 
theoretical work by Michel in 1964 [15] and that of other authors in the late 1960's 
[16, 17, 18], the results of such experiments have generally been analyzed in terms of a 
meson-exchange picture, and in 1980 the work of Desplanques, Donoghue, and Holstein 
(DDH) developed a comprehensive and general meson-exchange framework for the anal- 
ysis of such interactions in terms of seven parameters representing weak parity- violating 
meson-nucleon couplings [19]. The DDH interaction has become the standard setting by 
which hadronic and nuclear PV processes are now analyzed theoretically. 

It is important to observe, however, that the DDH framework is, at heart, a model 
based on a meson-exchange picture. Provided one is interested primarily in near-threshold 
phenomena, use of a model is unnecessary, and one can instead represent the PV nuclear 
interaction in a model-independent effective-field-theoretic fashion. The purpose of the 
present work is to formulate such a systematic, model-independent treatment of PV NN 
interactions. We feel that this is a timely goal, since such PV interactions are interesting 
not only in their own right but also as effects entering atomic PV measurements [21] 
as well as experiments that use parity violation in electromagnetic interactions to probe 
nucleon structure [23]. 
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In our reformulation of nuclear PV, we consider two versions of EFT, one in which the 
pions have been "integrated out" and the other including the pion as an explicit degree 
of freedom. In the pionless theory, the PV nuclear interaction is entirely short-ranged, 
and the most general potential depends at leading order on five independent operators 
parameterized by a set of five a priori unknown low-energy constants (LECs). When 
applied to low-energy (E cm < 50 MeV) two-nucleon PV observables — such as the neutron 
spin asymmetry in the capture reaction n + p — > d + 7 — it implies that there are five 
independent PV amplitudes, which may be determined by an appropriate set of mea- 
surements. We therefore recover previous results obtained without effective field theory 
by Danilov [24] and Desplanques and Missimer [25]. Making contact with these known 
results is an important motivation for us to consider this pionless EFT. Going beyond 
this, in next (non-vanishing) order in the EFT, there are several additional independent 
operators. By contrast, the DDH meson-exchange framework amounts to a model in 
which the short-range physics is codified into six independent operators. On one hand, 
the heavy-meson component of the DDH potential is a redundant representation of the 
leading-order EFT. On the other, it does not provide the most complete parameterization 
of the short-ranged PV NN force to subleading order, because it is based on a truncation 
of the QCD spectrum after inclusion of the lowest-lying octet of vector mesons. It may, 
therefore, not be entirely physically realistic, and we feel that a more general treatment 
using EFT is warranted. 

When we are interested in observables at higher energies, we need to account for pion 
propagation explicitly, simultaneously removing its effects from the contact interactions. 
Inclusion of explicit pions introduces a long-range component into the PV NN interaction, 
whose strength is set at the lowest order by the PV 7rNN Yukawa coupling, h\ NN . This 
long-range component, which is formally of lower-order than shorter-range interactions, is 
identical to the long-range, one-pion-exchange (OPE) component of the DDH potential. 
However, in addition, inclusion of pions leads to several new effects that do not arise 
explicitly in the DDH picture: 

• A medium-range, two-pion-exchange (TPE) component in the potential that arises 
at the same order as the leading short-range potential and that is also proportional 
to h\ NN . This medium-range component was considered some time ago in Ref. 
[18] but could not be systematically incorporated into the treatment of nuclear PV 
until the advent of EFT. As a result, such piece has not been previously included 
in the analysis of PV observables. We find that the two-pion terms introduce a 
qualitatively new aspect into the problem and speculate that their inclusion may 
modify the h\ NN sensitivity of various PV observables. 

• Next-to-next-to-leading-order (NNLO) PV irNN operators. In principle, there ex- 
ist several such operators that contribute to the PV NN interaction at the same 
order as the leading short-range potential. In practice, however, effects of all but 
one of the independent NNLO PV nNN operators can be absorbed via a suitable 
redefinition of the short-range operator coefficients and h\ NN . The coefficient of 
the remaining, independent NNLO operator - k]^ NN - must be determined from ex- 
periment. Additional terms are generated in the potential at 0(Q) by higher-order 
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corrections to the strong 71NN coupling (here, Q denotes a small momentum or pion 
mass). These terms have also not been included in previous treatments of the PV 
NN interaction. Their coefficients are fixed by either reparameterization invariance 
or measurements of other parity-conserving pion-nucleon observables. 

• A new electromagnetic operator. For PV observables involving photons, the explicit 
incorporation of pions requires inclusion of a PV NNir / ~f operator that is entirely 
absent from the DDH framework and whose strength is characterized by a constant 

In short, for the low-energy processes of interest here, the most general EFT treat- 
ment of PV observables depends in practice on eight a priori unknown constants when 
the pion is included as an explicit degree of freedom: five independent combinations of 
0(Q) short-range constants and those associated with the effects of the pion: h\ NN , C w , 
and the NNLO PV nNN coupling k\ a NN . In order to determine these PV low-energy 
constants (LECs), one therefore requires a minimum of five independent, low-energy ob- 
servables for the pionless EFT and eight for the EFT with dynamical pions. Given the 
theoretical ambiguities associated with interpreting many-body nuclear observables (see 
below), one would ideally attempt to determine the PV LECs from measurements in few- 
body systems. Indeed, the state of the art in few-body physics allows one perform ab 
initio computations of few-body observables [1], thereby making the few-body system a 
theoretically clean environment in which to study the effects of hadronic PV. At present, 
however, there exist only two measurements of few-body PV observables: A pp , the longi- 
tudinal analyzing power in polarized proton-proton scattering, and A 1 ^, the longitudinal 
analyzing power for pa scattering. In what follows, we outline a prospective program of 
additional measurements that would afford a complete determination of the PV LECs 
through 0(Q). 

Completion of this low-energy program would serve two additional purposes. First, it 
would provide hadron structure theorists with a set of benchmark numbers that are in 
principle calculable from first principles. This situation would be analogous to what one 
encounters in chiral perturbation theory for pseudoscalar mesons, where the experimental 
determination of the ten LECs appearing in the 0(Q A ) Lagrangian presents a challenge 
to hadron-structure theory. While many of the 0(Q 4 ) LECs are saturated by t-channel 
exchange of vector mesons, it is not clear a priori that the analogous PV NN constants are 
similarly saturated (as is assumed implicitly in the DDH model). Moreover, analysis of the 
PV NN LECs involves the interplay of weak and strong interactions in the strangeness- 
conserving sector. A similar situation occurs in AS = 1 hadronic weak interactions, and 
the interplay of strong and weak interactions in this case is both subtle and only partially 
understood, as evidenced, e.g., by the well-known AI =1/2 rule enigma. The additional 
information in the AS = sector provided by a well-defined set of experimental numbers 
would undoubtedly shed light on this fundamental problem. 

The information derived from the low-energy few-nucleon PV program could also pro- 
vide a starting point for a reanalysis of PV effects in many-body systems. Until now, one 
has attempted to use PV observables obtained from both few- and many-body systems 
in order to determine the seven PV meson-nucleon couplings entering the DDH potential, 
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and several inconsistencies have emerged. The most blatant is the vastly different value 
for h\ NN obtained from the PV 7-decays of 18 F and from the combination of the pp asym- 
metry and the Cesium anapole moment. Although combinations of coupling constants 
can be found that fit partial sets of experiments (see, e.g., Ref. [20]), it seems difficult 
to describe all experiments consistently with theory (see, e.g., Ref. [21] and references 
therein). The origin of this clash could be due to any one of a number of factors. Using 
the operator constraints derived from the few-body program as input into the nuclear 
analysis could help clarify the situation. It may be, for example, that the medium-range 
TPE potential or higher-order operators relevant only to nuclear PV processes play a more 
significant role in nuclei than implicitly assumed by the DDH framework. Alternatively, 
the treatment of the many-body system — such as the truncation of the model space in 
shell-model approaches to the Cesium anapole moment — may be the culprit. (For an 
example of the relevance of nucleon-nucleon correlations to parity violation in nuclei, see 
Ref. [22].) In any case, approaching the nuclear problem from a more systematic perspec- 
tive and drawing upon the results of few-body studies would undoubtedly represent an 
advance for the field. 

In the remainder of the paper, then, we describe in detail the EFT reformulation of 
nuclear PV and the corresponding program of study. In Section 2, we briefly review the 
conventional, DDH analysis and summarize the key differences with the EFT approach. 
In particular, we write down the various components of the PV EFT potential here, rele- 
gating its derivation to subsequent sections. In Section 3, we outline the phenomenology 
of the low-energy few-body PV program, providing illustrative relationships between vari- 
ous observables and the five relevant, independent combinations of short-range LECs. We 
emphasize that the analysis presented in Section 3 is intended to demonstrate how one 
would go about carrying out the few-body program rather than to give precise numerical 
formulas. Obtaining the latter will require more sophisticated few-body calculations than 
we are able to undertake here. Section 4 contains the derivation of the PV potential in 
the EFT without explicit pions. We then extend the framework to include pions explicitly 
in Section 5. In Section 6 we discuss the relationship between the PV LECs and the PV 
meson-nucleon couplings entering the DDH framework, and illustrate how this relation- 
ship depends on one's truncation of the QCD spectrum. Section 7 contains some final 
observations. Various details pertaining to the calculations contained in the text appear 
in the Appendices. 



The essential idea behind the conventional DDH framework relies on the fairly success- 
ful representation of the parity-conserving NN interaction in terms of a single meson- 
exchange approach. Of course, this requires the use of strong-interaction couplings of the 
lightest vector (p, a;) and pseudoscalar (tt) mesons(M), 
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Figure 1: Parity- violating NN potential generated by meson exchange. 



whose values are reasonably well determined. The DDH approach to the parity-violating 
weak interaction utilizes a similar meson-exchange picture, but now with one strong and 
one weak vertex — cf. Fig. 1. 

We require then a parity- violating NNM Hamiltonian in analogy to Eq. (1). The pro- 
cess is simplified somewhat by Barton's theorem, which requires that in the CP-conserving 
limit, which we employ, exchange of neutral pseudoscalars is forbidden [26]. From general 
arguments, the effective Hamiltonian with fewest derivatives must take the form 
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wk = ^N(t x n) 3 N - N (h° p r ■ (f + h\p% + - r ■ ff)^ W N 
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We see that there exist, in this model, seven unknown weak couplings h\ NNl h p , ... How- 
ever, quark model calculations suggest that h' p l is quite small [27], so this term is usually 
omitted, leaving parity- violating observables described in terms of just six constants. 
DDH attempted to evaluate such PV couplings using basic quark-model and symmetry 
techniques, but they encountered significant theoretical uncertainties. For this reason 
their results were presented in terms of an allowable range for each, accompanied by a 
"best value" representing their best guess for each coupling. These ranges and best values 
are listed in Table 1, together with predictions generated by subsequent groups [28, 29]. 

Before making contact with experimental results, however, it is necessary to convert 
the NNM couplings generated above into a parity-violating NN potential. Inserting the 
strong and weak couplings, defined above, into the meson-exchange diagrams shown in 
Fig.l and taking the Fourier transform, one finds the DDH parity- violating NN potential 
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Table 1: Weak NNM couplings as calculated in Refs. [19, 28, 29]. All numbers are 
quoted in units of the "sum rule" value g w = 3.8 ■ 10~ 8 . 
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where pi = —iVi, Vj denoting the gradient with respect to the co-ordinate Xj of the i-th 
nucleon, r = \x± — x 2 \ is the separation between the two nucleons, 



Wi(r) 



exp(— rriir) 



(4) 



is the usual Yukawa form, and the strong ttNN coupling g w NN has been expressed in 
terms of the axial-current coupling qa using the Goldberger-Treiman relation: g^NN = 
gAm^/Fir, with F n = 92 A MeV being the pion decay constant. 

Nearly all experimental results involving nuclear parity violation have been analyzed 
using VpE, H for the past twenty-some years. At present, however, there appear to exist 
discrepancies between the values extracted for the various DDH couplings from experi- 
ment. In particular, the values of h\ NN and h° p extracted from pp scattering and the 7 
decay of 18 F do not appear to agree with the corresponding values implied by the anapole 
moment of 133 Cs measured in atomic parity violation [30]. 

These inconsistencies suggest that the DDH framework may not, after all, adequately 
characterize the PV NN interaction and provides motivation for our reformulation using 
EFT. The idea of using EFT methods in order to study parity-violating hadronic AS = 
interactions is not new [31]. Recently, a flurry of activity (see, for example, Refs. 
[32, 33, 34, 35, 36, 37, 38]) has centered on PV processes involving a single nucleon, such 
as 

ep — > ep, 7p — > 7p, TP ~ etc. 

There has also been work on the iViV system, with pion exchange treated perturbatively 
[39, 40] or non-perturbatively [41]. However, a comprehensive analysis has yet to take 
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place, and this omission is rectified in the study described below, wherein we generate 
a systematic framework within which to address PV NN reactions. We utilize the so- 
called Weinberg formulation [3], wherein the pion, when included explicitly, is treated 
fully quantum mechanically while shorter-distance phenomena — as would be produced 
by the exchange of heavier mesons such as p, u, etc. — are represented in terms of simple 
four-nucleon contact terms. The justification for the non-perturbative treatment of (parts 
of) pion exchange has been discussed in a recent paper [42]. 

Although a fully self-consistent procedure would involve use of EFT to compute both 
the PV operators and few-body wavefunctions, equally accurate results can be obtained by 
drawing upon state-of-the art wave functions obtained from a phenomenological, strong- 
interaction A^ A" potential, including PV effects perturbatively, and using EFT to system- 
atically organize the relevant PV operators. Such a "hybrid" approach has been followed 
with some success in other contexts [2] and we adopt it here. In so doing, we truncate our 
analysis of the PV operators at order Q/A x , where Q is a small momentum characteristic 
of the low-energy PV process and A x = 4irF n ~ 1 GeV is the scale of chiral symmetry 
breaking. Since realistic wave functions obtained from a phenomenological potential effec- 
tively include strong-interaction contributions to all orders in Q/A x , the hybrid approach 
introduces some inconsistency at higher orders in Q/A x . For the low-energy processes 
of interest here (E cm < 50 MeV), however, we do not expect the impact of these higher- 
order problems to be significant. We would not, however, attempt to apply our analysis to 
higher-energy processes (e.g., the TRIUMF 221 MeV pp experiment [43]) where inclusion 
of higher-order PV operators would be necessary. 

With these considerations in mind, it is useful to compare V^vk with the leading-order 
PV NN EFT potential. In the pionless theory, this potential is entirely short ranged 
(SR) and has co-ordinate space form 



A 3 

x 



Ci + (C 2 + C 4 ) (p^p) 3 + C 3 n ■ r 2 + X ah C h Tlrl 

(<ti - <r 2 ) • {-iV,f m (r)} 
+ [Ci + (C 2 + C< 4 ) (^ip) 3 + <Vi • r 2 + l ab C 5 T?r> 

i (<?i x <j 2 ) • [-iV, f m {r)\ 
'Ti - r 2 



+ (C 2 - C 4 ) 



ffi + o 2 ) ■ {-iV, f m {r)} 



+C & ie ah ^T b 2 («?! + a 2 ) ■ HV, f m (r)}} 
where the subscript "1" in the potential denotes the chiral index of the operators 1 , 



(5) 




(6) 



and f m (r) is a function that 

1 Roughly speaking, the chiral index corresponds to the order of a given operator in the Q/A x expan- 
sion. A precise definition is given in Sec. 5 below. 
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i) is strongly peaked, with width ~ 1/m about r = 0, and 

ii) approaches 5^(f) in the zero-width (m — > oo) limit. 
A convenient form, for example, is the Yukawa-like function 

2 

Tfl 

fm(r) = — exp (-mr). (7) 

Here m is a mass chosen to reproduce the appropriate short-range effects (m ~ A x in the 
pionful theory, but m ~ m n in the pionless theory). Note that, since the terms containing 
C2 and C4 are identical, nominally contains ten independent operators. As we show 
below, however, only five combinations of these operators are relevant at low-energies. 

For the purpose of carrying out actual calculations, one could just as easily use the 
momentum-space form of V^ P ^ R , thereby avoiding the use of f m {r) altogether. Neverthe- 
less, the form of Eq. (5) is useful when comparing with the DDH potential. For example, 
we observe that the same set of spin-space and isospin structures appear in both V^ P ^ R 
and the vector-meson exchange terms in V^^, though the relationship between the vari- 
ous coefficients in V P g R is more general. In particular, the DDH model is tantamount to 
taking m ~ m p , and assuming 

Ci = c£ Xh " ^' 

— — — — — — 1 4- (Q\ 

<^3 ^4 ^5 

assumptions which may not be physically realistic. In Section 6, we give illustrative 
mechanisms which may lead to a breakdown of these assumptions. 

When pions are included explicitly, one obtains in addition the same long-range (LR) 
component induced by OPE as in I^dhj 



3 



P1-P2 f X 

wJr) 



2m 



N 



(10) 



where w n (r) is given by Eq. (4). Note that, as we will explain in Sect. 5, V^ LR is 
two orders lower than V± g R — in contrast to the strong potential where the short- and 
long-range components first arise formally at the same order. (Even though Eq. (10) has 
the same form as a term in Eq. (5), it has no suppression by powers of A x or other heavy 
scales. Therefore, it appears at lower order.) 

Furthermore, two new types of contributions to the potential arise at the same order 
as Eq. (5): (a) a long-range component stemming from higher-order nNN operators, 
and (b) a medium-range (MR), two-pion- exchange (TPE) contribution, Vf Y MR . At 0(Q), 
the TPE potential is proportional to h\ NN and involves two terms having the same spin- 
isospin structure as the terms in proportional to C* 2 and Cq but having a more 

complicated spatial dependence. In momentum space, 

V&M = ~ {C?(q)^^i (*i x a 2 ) ■ q + Cf (q)ie ab3 r^ (a, + a 2 ) ■ q] , (11) 
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where the functions C$*(q) and C^i^q), defined below in Eq. (121), are determined by the 
leading-order tiNN couplings. Again, it is more convenient to compute matrix elements 
of V\ P ^[ R using the momentum-space form, and we defer a detailed discussion of the 
latter until Section 5 below. We emphasize, however, the presence of V 1 p 'm R introduces a 
qualitatively new element into the treatment of nuclear PV with pions not present in the 
DDH framework. 

The NNLO long-range contribution to the potential generated by the new PV irNN 
operator is 



^ PV lr = -*f^f (^P) 3 {e abc ^ 2 e {V?, [V>V^(r)]} + (1 ~ 2)} + • ■ • , (12) 



where V r is the gradient with respect to the relative co-ordinate x\ — x 2 and where the 
• • • denote long-range, NNLO contributions proportional to h\ NN that are generated by 
NNLO effects at the strong irNN vertex (see Appendix B). 

As we discuss in Section 3, a complete program of low-energy PV measurements in- 
cludes photo-reactions. In the DDH framework, PV electromagnetic (EM) matrix el- 
ements receive two classes of contributions: (a) those involving the standard, parity- 
conserving EM operators in combination with parity-mixing in the nuclear states, and 
(b) PV two-body EM operators derived from the amplitudes of Fig. 13. Explicit expres- 
sions for these operators in the DDH framework can be found in Ref. [21]. In the case of 
EFT, the two-body PV EM operators associated with heavy-meson exchange in DDH are 
replaced by operators obtained by gauging the derivatives in V 1 F g R as well as by explicit 
photon insertions on external legs. The two-body operators associated with LR are 
identical to those appearing in DDH, while the PV currents associated with V^\ AR and 
Vi, lr are obtained by gauging the derivatives appearing in the potential and by inserting 
the photon on all charged-particle lines in the corresponding Feynman diagrams 2 . The 
foregoing two-body currents introduce no new unknown constants beyond those already 
appearing in the potential. However, an additional, independent pion-exchange two-body 
operator also appears at the same order as the short-range PV two-body currents: 

J{2i,32,$ = ^gffi^ e-** r+r 2 - a, x q a 2 ■ r h n (r) + (1 <- 2), (13) 

where 

, . . exp (— m^r) ( 1 \ 

K(r) = — ^ 1 + , 14 

m n r V m w r J 

and CV is an additional LEC parameterizing the leading PV NNir / ~f interaction. Any 
photoreaction sensitive to the short-range PV potential will also depend on C n when 
pions are included explicitly. 

Through 0(Q), then, the phenomenology of nuclear PV depends on five unknown 
constants in the pionless theory and eight when pions are included explicitly {h^N, 
k\ a NNl and C w in addition to the contact interactions). As we discuss below, an initial 

2 The derivation of the medium-range two-body operators involves an enormously detailed computa- 
tion, which we defer to a later publication. 
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low-energy program will afford a determination of the five constants in the pionless theory. 
Additional low-energy measurements in few-body systems would provide a test of the 
self-consistency of the EFT at this order. Any discrepancies could indicate the need to 
including pions as explicit degrees of freedom, thereby necessitating the completion of 
additional measurements in order to determine the pion contributions to 0(Q). As we 
discuss below, there exists a sufficient number of prospective measurements that could 
be used for this purpose. Given the challenging nature of the experiments, a sensible 
strategy would be to first test for the self-consistency of the pionless EFT with a smaller 
set of measurements and then to complete the additional measurements needed for EFT 
with pions if necessary. 

3 Parity Violation in Few-Body Systems 

There exist numerous low-energy experiments that have attempted to explore hadronic 
parity violation. Some, like the photon asymmetry in the decay of a polarized isomeric 
state of 180 Hf, 

A 1 = -(1.66 ±0.18) x l(r 2 [44], (15) 
or the asymmetry in longitudinally-polarized neutron scattering on 139 La, 

A z = (9.55 ± 0.35) x 10~ 2 [45], (16) 

involve F-P shell nuclei wherein the effects of hadronic parity violation are large and 
clearly observed, but where the difficulty of performing a reliable wave function calculation 
precludes a definitive interpretation. For this reason, it is traditional to restrict one's 
attention to S-D shell or lighter nuclei. Here too, there exist a number of experiments, 
such as the asymmetry in the decay of the polarized first excited state of 19 F, 

Ayd ,110keV) = -(8.5 ±2.6) x 10~ 5 [46] 

= -(6.8 ± 1.8) x 10~ 5 [47], (17) 

wherein a clear parity- violating signal is observed, or those such as the circular polarization 

(24 ±24) x HT 4 [48] 

(3± 16) x 10~ 4 [49], (18) 

(-7 ±20) x 10~ 4 [50] 
(3 ±6) x 10~ 4 [51] 
(—10 ± 18) x 10" 4 [52] 
(2 ±6) x 10" 4 [53], (19) 



in the decay of excited levels of Ne, 

P 7 (- ,2.789 MeV) 
2 



or of 18 F, 

P 7 (0~, 1.081 MeV) 
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wherein a nonzero signal has not been seen, but where the precision of the experiment 
is high enough that a significant limit can be placed on the underlying parity-violating 
mechanism. 

The reason that a 1CT 4 experiment can reveal information about an effect which is on 
the surface at the level 



where pp ~ 270 MeV is the Fermi momentum, is that the nucleus can act as an PV- 
amplifier. This occurs when there exist a pair of close-by levels having the same spin 
but opposite parity, (J 1 * 1 ). In this case the parity mixing expected from lowest-order 
perturbation theory, \ip±) ~ | J^) ± e\ J T ), can become anomalously large due to the 
smallness of the energy difference Ej+ — Ej- in the mixing parameter 



Indeed, when compared with a typical level splitting of ~ 1 MeV, the energy differences 
exploited in 19 F (AE = 110 keV), 21 Ne (AE = 5.7 keV), and 18 F (AE = 39 keV) lead 
to expected enhancements at the level of 10, 100, and 25 respectively. However, when 
interpreted in terms of the best existing nuclear shell-model wave functions, there exists 
a serious discrepancy between the values of the AI=1 pion coupling required in order to 
understand the 19 F or 21 Ne experiments and the upper limit allowed by the 18 F result. 

Such matters have been extensively reviewed by previous authors [54, 55, 56], and 
we do not intend to revisit these issues here. Instead we suggest that at the present 
time any detailed attempt to understand the parity-violating NN interaction must focus 
on experiments involving only the very lightest — NN, Nd, Na — systems, wherein our 
ability to calculate the effects of a given theoretical picture are under much better control. 
As we demonstrate below, there exist a sufficient number of such experiments, either in 
progress or planned, in order to accomplish this task for either the pionless EFT or the 
EFT with dynamical pions. Once a reliable set of low-energy constants are in hand, 
as obtained from such very-light systems, theoretical work can proceed on at least two 
additional fronts: 

i) experimental results from the heavier nuclear systems — involving P, S-D, and F-P 
shells and higher levels — can be revisited and any discrepancies hopefully resolved 
with the confidence that the weak low-energy constants are correct, and 

ii) one can attempt to evaluate the size of the phenomenological weak constants starting 
from the fundamental quark-quark weak interaction in the Standard Model. 

This scheme mirrors the approach that has proven highly successful in chiral per- 
turbation theory (ChPT) [57], wherein phenomenological constants are extracted purely 
from experimental results, using no theoretical prejudices other than the basic (broken) 
chiral symmetry of QCD. In the meson sector [58], the phenomenologically-determined 
counter-terms L\, L 2 , L w have already become the focus of various theoretical pro- 
grams attempting to predict their size from fundamental theory. Note that our approach 
to nuclear parity violation is similar in spirit to the one advocated in a prescient 1978 



G^m^ x {pf/^n) ~ 10 




(20) 
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paper by Desplanques and Missimer [25] that builds on ideas put forward by Danilov [24] . 
In subsequent work, this approach was superseded by the use of the DDH potential. In 
our study, then, we are in a sense recasting the ideas of Refs. [25, 24] in the modern and 
theoretically systematic framework of EFT. 



3.1 Amplitudes 

We now consider the first part of the program — elucidation of the basic weak couplings. 
We argue that, provided one is working in a region of sufficiently low-energy, the parity- 
violating NN interaction can be described in terms of just five real numbers, which 
characterize S-P wave mixing in the spin singlet and triplet channels. The arguments in 
this section borrow heavily from the work of Danilov [24] and Desplanques and Missimer 
[25]. The following sections will show how to interpret this phenomenology within EFT. 

For simplicity we begin with a parity-conserving system of two nucleons. Then the NN 
scattering matrix can be expressed purely in terms of S-wave scattering at low energies 
and has the phenomenological form [24] 

M(k f , = (k f \f\ki) = m t (A;)Pi + m s (k)P , (21) 

where 

Pi = -(3 + <7i • <r 2 ), P = -(l-Bi-B 2 ) 

are spin-triplet, -singlet projection operators. All other partial waves can be neglected. 
We can determine the form of the functions rriiik) by using the stricture of unitarity, 

2Imf = T f T. (22) 

In the S-wave sector this becomes 

Immi(fe) = k\ mi {k)\ 2 , (23) 
whose solution is of the familiar form 

mi (k) = le iSiW sin Si{k). (24) 

Since at zero energy 

limmjf&O = — a,i (25) 

where Oj is the scattering length, it is clear that unitarity can be enforced by the simple 
modification 

mik) = — (26) 
1 + ikdi 

which is the lowest-order effective-range result. The scattering cross section is found via 

% = M = TTW' (27) 
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so that at the lowest energy we have the familiar form 

da s t , ,o 



lim 



a s,t 



k-^0 dfl 

The associated scattering wave functions are given by 

ik\r—r'\ 



(28) 



4 +) (f) = 



ik-f 



47T 



jdh' B -^- V{r')^\r') 



X 



„ikr 



e lk - r + M(-iV \k)- 



X, 



(29) 



where \ is the spin wave function. In the simple Born approximation, then, we can 
represent the wave function in terms of an effective local potential 



v (0)PC 



(f) 



An 

m N 



(a t Pi + a s P )<f 3 (r), 



(30) 



as can be confirmed by substitution into Eq. (29). 

Parity mixing can be introduced into this simple representation, as done by Danilov 
[24], via generalization of the scattering amplitude to include parity- violating structures. 
Up to laboratory momenta of 140 MeV or so, we can omit all but S- and P-wave mixing, 
in which case there exist five independent such amplitudes: 

i) d t (k), representing 3 S\ — l P\ mixing; 

ii) d° s ,1,2 {k), representing 1 S — 3 P mixing with AI = 0, 1, 2 respectively; and 

iii) Ct(k), representing 3 Si — 3 Pi mixing. 

and, after a little thought, it becomes clear that the low-energy scattering matrix in the 
presence of parity violation can be generalized to 

M(k f ,ki) = m t (k)P 1 + m s (k)P 

+ [(d° s (k)Q 1 + d\{k)Q 1+ + d 2 s (k)Q 2 ) (h • («?! - a 2 )P 1 + pSj ■ (<?i - a 2 )) 

+d t (k) fa • (^ - a 2 )P + P k f • (*! - a 2 )) 
+c t {k)Q 1 _{& 1 + a 2 ) ■ (hP, + P-ykf) , (31) 

where we have introduced the isovector and isotensor operators 



-{Tl-T 2 ) z , 



Ql+ = ^(Tl + T 2 ) z , 



2V6 



(3t 1z t 2z - fi • f 2 ), (32) 



and isospin projection operators 



Qo = J(l -Ti -T 2 ), Qi = + -T 2 ). 



(33) 
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Each of the new pieces is indeed odd under spatial inversion [c?j — > <jj and kf,ki — >■ 
— fc/, — fcj] and even under time reversal — > — <7j and fcj • (<ti — <J 2 )Pj <-> -P?^/ • (<?i — 0*2)]- 
Now consider what constraints can be placed on the forms di(k), Ci(k). The requirement 
of unitarity reads 

Im di(k) = k[m*(k)di(k) + d*(k)m p (k)}, (34) 

where rrii(k) , m p (k) are the scattering amplitudes in the S-, P-wave channels connected 
by di(k). Eq. (34) is satisfied by the solution 

di(k) = |^(A;)|expi[^(A;) + 5 p (k)}, (35) 

i.e., the phase of the transition amplitude is simply the sum of the strong interaction 
phase shifts in the incoming and outgoing channels. 

Danilov [24] suggested that, on account of the short-range of the weak interaction, the 
energy dependence of the weak couplings di(k) should be primarily determined, up to say 
50 MeV or so solely by the strong interaction dynamics. Since at very low energy the 
P-wave scattering can be neglected, he suggested the use of the forms 

c t (k) = p t m t (k), d t (k) = X t m t (k), dl(k) = Km s (k), (36) 

which provide the parity-mixing amplitudes in terms of the five phenomenological con- 
stants: p t , At, A* . 

We can understand the motivation behind Danilov's assertion by writing down the 
simplest phenomenological form for a weak low-energy parity-violating NN potential. To 
do so, one may start with the momentum-space form of V^g R given in Eq. (5): 



A 3 

x 



(Ci + C 3 )Qi + (Ci - 3G0Q o + (C 2 + C 4 )Q 1+ - d^C 5 Q 2 



+ 



(Ci + C 3 )Q 1 + (d - 3C 3 )g + (C 2 + C 4 )Q 1+ - J^C 5 Q 2 



i (<7i x a 2 ) ■ q 

+ [C 2 - C 4 ] Qi_ (ai + a 2 ) ■ p + C 6 i (n x r 2 ) z (^ + a 2 ) ■ q\, (37) 



where p = [(pi - p 2 ) + (pi - p' 2 )]/2 and g = [(pi - p 2 ) - (p[ - p' 2 )]/2. 

The change in the wave function generated by Vi \ R (r) is understood to involve the 
full strong-interaction Green's function and wave functions 

^(+)(f) = J d 3 r'G k (f^)V^l R (^ +) ^) (38) 

and the connection between the weak PV potential V^ R (r) and the scattering matrix 
Eq. (31) can be found via 



di(k) 



mK (<tf ( -fete S(+, > + (^-'IV^kI^ (+) )). 



4:71 



(39) 
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Now, if we are at very low energy, we may use the plane-wave approximation for the P 

wave, 

rf'W^JiW, (40) 
and we can approximate the S wave by its asymptotic form 

+ (r) ~ -L e *« sin(Ax + ^(jfe)) ^ -L e i5 ^ sin Si(k). (41) 

KT KT 

where we have used the experimental fact that \Si\ ~ \kai\ » k/m where 1/m is max- 
imum range set by the integration. Then, we can imagine calculating a generic parity- 
violating amplitude di(k) via Eq. (39): 



4.77- roc f) 1 

d t {k) ~ T Qj o drr^ikr^JUr)}^ 6 ^ sin 5 t (k) 
1 



Aj^e i(5i sin^j = \rrii(k) (42) 

with "Cj" symbolically indicating the appropriate combination of PV constants appearing 
in V\ pi g R and 

4tt roo ■ v djrnir) 4tt /■«> 2 rf/ m (r) 
Ai ~ -pCj / drrji{kr)— — « — C { \ drr — — , (43) 
k Jo dr 3 Jo dr 

which is the basic form advocated by Danilov 3 . An analogous relationship holds for p t . 

At low energy then it seems prudent to explicitly include the appropriate S-wave scat- 
tering length in expressing the effective weak potential, and we can define 

\imm s +(k) = -a S)t , lim c t (k), d s (k), d t (k) = -p t a t , -A*a s , -X t a t (44) 

fe^O k^O 

As emphasized above, the real numbers p t , A*, A t — which can in turn be related to the 
effective parameters Cj — completely characterize the low-energy parity-violating interac- 
tion and can be determined experimentally, as we shall discuss below. Alternatively, we 
use an isospin decomposition 

K p = ^s-^l 

K n = \° s - \l + (45) 

to write things in terms of the appropriate NN quantities. Now, the explicit connection 
between the S-matrix elements \, p t and the weak interaction parameters Cj,(7j in our 



3 Note that this argument is not quite correct quantitatively. Indeed, since 1/m is much smaller than 
the range of the NN interaction (in the pionful theory), we should really use interior forms of the wave 
functions. However, when this is done, the same qualitative result is found, but the simple relationship 
in Eq. (43) is somewhat modified. 
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effective Lagrangian must be done carefully using the Eq. (39) and the best possible NN 
wave functions. This work is underway, but is not yet completed [59]. In the meantime, 
we may obtain an indication of the connection by using the following simple arguments: 
When we restrict ourselves to a model-space containing only the low-energy S,P am- 
plitudes noted above, then several of the operators in Eqs. (5,37) become redundant. For 
example, the d t amplitude involves aT = 0->T = transition, so only the terms propor- 
tional to Qo contribute. In this case, the spin-space operators (<7i — a 2 ) -p and x a 2 ) • Q 
yield identical matrix elements up to an overall constant of proportionality. This feature 
can be seen by considering the co-ordinate space potential, which contains the function 
f m (r) times derivatives acting on the initial and final states. In the short range limit and 
in the absence of the NN repulsive core, both the P-wave and first derivative of the S-wave 
vanish at the origin, whereas the product of the S-wave and first derivative of the P-wave 
are non-zero. Thus, only the components of q and p that yield derivatives of the P-wave 
at the origin contribute, leading to identical matrix elements of these two operators (up to 
an overall phase). Of course, corrections to this statement occur when f m (r) is smeared 
out over some short range ~ 1/m. Since 1/m « 1/typical momentum ~ a where a is the 
scattering length, at low energy such corrections are higher-order in our power counting, 
going as K 2 /m 2 , where K ~ \J M(E + V) with V ~ 50 MeV representing some average 
depth of the NN potential characterizing the interior region. Similarly, the operators 
(<7i — (72)2 and (<7i x a 2 ) z each transform a spin triplet into a spin singlet state, and vice 
versa. Hence, one may absorb the effect of the term proportional to {C\ + C 3 ) into the 
corresponding term proportional to (C\ + C 3 ) by a suitable redefinition of the constants. 
Related arguments allow one to absorb the remaining terms proportional to the Ci - as 
well as the term containing Cq - into the terms involving [C\ — 3Cs)P\, (C 2 + Ci)Q\ + , 
(C 2 — C^)Qi-, and C^Q 2 for a net total of five independent operators, which in turn gen- 
erate the five independent low-energy PV amplitudes A°, A^, A^, A t , p t . In the zero-range 
limit, then, we have 

A t oc (d - 3C 3 ) - (Cx - 3C 3 ) 
A° oc (Ci + C 3 ) + (C 1 ! + C 3 ) 

\] oc (C 2 + C 4 ) + (C 2 + C4) (46) 
A^ oc -^8/3(C 5 + C 5 ) 

Pt oc l -(C 2 -C A )-C % . 

However, going away from strict threshold values and the use of more realistic wave 
functions will modify these expectations somewhat, as illustrated by a simple, didactic 
discussion in Appendix E. We emphasize, however, that what is needed at the present 
time is a purely empirical evaluation in terms of five independent and precise experiments, 
and that is what we shall discuss next. 

3.2 Relation to observables 

The next step of the program — contact between this effective parity- violating interaction 
and experimental observables — was initiated by Desplanques and Missimer [25]. Before 
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quoting these results, we sketch the manner by which such a confrontation is performed. In 
doing so, we emphasize that the following analysis does not rely on definitive computations 
employing state-of-the art few-body wave functions — carrying out such calculations goes 
beyond the scope of the present study. Indeed, obtaining precise values for the Aj and 
p t will require a concerted effort on the part of both experiment and few-body nuclear 
theory. What we provide below is intended, rather, to serve as a qualitative roadmap for 
such a program, setting the context for what we hope will be future experimental and 
theoretical work. 

For simplicity, we begin with an illustrative example of nn scattering, for which the 
Pauli principle demands that the initial state at low energy must be purely 1 Sq. One 
can imagine longitudinally polarizing a neutron of momentum p and measuring the total 
scattering cross section from an unpolarized target. Since a ■ p is odd under spatial 
inversion, the cross section can depend on helicity only if parity is violated, and via trace 
techniques the helicity-correlated cross section can easily be found. Using 

M(k f , h) = m s (k)P + d n s n [h • (<?! - a 2 )P + P k f • (<?! - «7 2 )] (47) 

we determine 

o~± — J <m f TrM(k f , h)\(l ± <?! • k)M\k f , h) 

= 4vr|m s (A;)| 2 TrPo + 8vrRe m*(k)<% n (k)Tr P (<?i - p 2 ) ■ k(l ± <?i • k) + • • • 

= 47r|m s (A;)| 2 ±167rRem:(A;X ri (A;) + --- (48) 

Defining the asymmetry via the sum and difference of such helicity cross sections and 
neglecting the tiny P-wave scattering, we have then 

a + -a_ 4kRe[m* s (k)d™ n (k)} A1 xnn , An . 

A L = — = ' , . ' J ~ 4/cA"". (49) 

o"+ + cr_ \m s (k)\ z 

Thus the helicity-correlated nn asymmetry provides a direct measure of the parity- 
violating parameter A™ n . Note that, since the total cross section is involved, some in- 
vestigators have opted to utilize the optical theorem via [60, 61, 62] 

which, using our unitarized forms, is completely equivalent to Eq. (49). 

Of course, nn scattering is currently just a gedanken experiment, and we have dis- 
cussed it merely as a warm-up to the real problem: pp scattering, which introduces the 
complications associated with the Coulomb interaction. In spite of this complication, the 
calculation proceeds quite in parallel to the discussion above with obvious modifications. 
We find 

^^yyii^. (51) 

a + + a- \m s (k)\ 2 

In the next section we show how this can be obtained straightforwardly within an EFT 
approach. 
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On the experimental side, such asymmetries have been measured both at low energies 
(13.6 and 45 MeV) as well as at higher energies (221 and 800 MeV). It is only the low- 
energy results 4 

Af (13.6 MeV) = -(0.93 ± 0.20 ± 0.05) x 10~ 7 [63], 
Af(45MeV) = -(1.50 ± 0.22) x 10" 7 [64], (53) 

that are appropriate for our analysis, and from these results we can extract the experi- 
mental number for the singlet mixing parameter as 

(Af )exp t = _ A L (45MeV) = _ (4 Q ± Q g) x 1Q _ 8 fm (M) 

where Ak 0.88m7v- Note that this Eq. (54) is consistent with that of Desplanques and 
Missimer[25] 

(Af )exp t = _ A L (45MeV) = _ (4 A ± Q g) x 1Q _ 8 ^ (55) 
0.82771^ 

In a corresponding fashion, as described by Ref. [25], contact can be made between 
other low-energy observables and the effective parity-violating interaction. Clearly we 
require five independent experiments in order to identify the five independent S-P mixing 
amplitudes. As emphasized above, we consider only PV experiments on systems with 
A = 4 or lower, in order that nuclear-model dependence be minimized. We utilize here 
the results of Desplanques and Missimer [25], but these forms should certainly be updated 
using state-of-the art few-body computations. There exist many such possible experiments 
and we suggest the use of 

i) low-energy pp scattering, for which 

pp(13.6 MeV) : Af = -0A8X p s p m N 
pp(45 MeV) : Af = -0.82Af m N (56) 

ii) low-energy pa scattering, for which 

pa(46 MeV) : A p L a = [-0.48(Af + ^Af ) - 1.07(p t + \\ t )}m N (57) 

iii) threshold np radiative capture for which there exist two independent observables 

circular polarization : P 7 = (0.63A t — 0.16A^ n )mjv 
photon asymmetry : A 1 = —0.107p t rriN (58) 

4 Note that the 13.6 MeV Bonn measurement is fully consistent with the earlier but less precise number 

Af (15 MeV) = -(1.7 ± 0.8) x 1(T 7 [65] (52) 

determined at LANL. 
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iv) neutron spin rotation from 4 He 



drh na 1 1 

-|- = [1.2(Ar + g AD - 2.68(p t - -\ t )]M n rad/m (59) 

Inverting these results, we can determine the five S-P mixing amplitudes via 

m N \PP = -1.22 (45 MeV) 

TfiNPt = — 9.35 A 1 (np — > dj) (60) 
mjvAf = 1.6 (45 MeV) -3.7^(46 MeV) + 37 A 7 (np^d 7 ) 
-2P 7 (np -> d 7 ) 

m^Aj = 0.4 (45 MeV) - 0.7 A p L a (46 MeV) + 7A 7 (np -> d 7 ) 
+ P 1 (np — > d 7 ) 

mjv A" n = 0.83 33.3 A^np -> d 7 ) - 0.69 (45 MeV) 

+ 1.18 (46 MeV) - 1.08P 7 (np -> d 7 ) (61) 

At the present time only two of these numbers are known definitively — the longitudinal 
asymmetry in pp, Eq. (53), and in pa scattering, 

(46 MeV) = -(3.3 ± 0.9) x 10 -7 [66]. (62) 

However, efforts are underway to measure the photon asymmetry in radiative np capture 
at LANSCE [67] as well as the neutron spin rotation on 4 He at NIST [68]. These measure- 
ments are also proposed at the neutron beamline at the Spallation Neutron Source (SNS) 
under construction at Oak Ridge National Laboratory. An additional, new measurement 
of the circular polarization in np radiative capture would complete the above program, 
although this is very challenging because of the difficulty of measuring the photon helicity. 
Alternatively, one could consider the inverse reaction — the asymmetry in 7 d — > np — and 
this is being considered at Athens [69] and at HIGS at Duke [70]. 

To the extent that one can neglect inclusion of the 7r as an explicit degree of freedom, 
one could use this program of measurements to perform a complete determination of the 
five independent combinations of 0(Q), PV LECs. Nonetheless, in order to be confident 
in the results of such a series of measurements, it is useful to note that other light systems 
can and should also be used as a check of the consistency of the extraction. There are 
various possibilities in this regard, including 

i) pd scattering 

A p L d (15 MeV) = (-0.21^ - 0.07Af - 0.13A t - 0.04Af )m N (63) 

ii) radiative nd capture 

A 7 = {\A2pt + 0.59A T ; n + 1.18A t + 0.51Af )m N (64) 
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iii) neutron spin rotation on H 




= (1.26/0* - 0.63A t + 1.8A T ; p + 0.45Af + 0A5\™)m N rad/m (65) 



Note that possible follow-ups of the LANSCE and NIST experiments include the last 
three processes [71]. 

We emphasize that the above results have been derived under the assumption that the 
spin-conserving interaction p t is short ranged - an assumption applicable at energies well 
below the pion mass. On the other hand, for the 46 MeV pa measurement, the proton 
momentum is well above m n , so integrating out the pion may not be justified. In this 
case, inclusion of the pion will lead to modification of the above formulas, introducing a 
dependence on h\ NN) k\ a NNl and C n . Thus, a total of eight low-energy few-body mea- 
surements would be needed to determine the relevant set of low-energy constants. In the 
foregoing discussion, we have identified eight few-body observables that could be used for 
this purpose. Additional possibilities include the PV asymmetry in near-threshold pion 
photo- or electro-production [37, 38, 73] or deuteron photodisintegration[74]. At present, 
we are unable to write down the complete dependence of the few-body PV observables on 
h\ NN) k\ a NN) and CV, since only the effects of the LO OPE PV potential (and associated 
two-body currents) have been included in previous few-body computations. Obtaining 
such expressions is a task requiring future theoretical effort. In any case, it is evident 
from our discussion that there exists ample motivation for several new few-body PV ex- 
periments and that a complete determination of the relevant PV low-energy constants is 
certainly a feasible prospect. 



4 EFT without explicit pions 

Although the foregoing analysis relied on traditional scattering theory, it is entirely equiv- 
alent to an EFT approach. In the following two sections, we present this EFT treatment 
in greater detail, considering first only processes where the momenta p of all external 
particles are much smaller than the pion mass. In this regime, the detailed dynamics 
underlying the NN interaction cannot be resolved, and interactions are represented by 
simple delta-function potentials. As with any EFT, this approximation is justified by 
a separation of scales. In this case, one has scales set by the NN scattering lengths 
— a s ~ —20 fm, a t ~ 5 fm — that are both much larger than the ~ l/m n range of the 
pion-exchange component of the NN strong interaction [8, 9]. Because of this separation 
of scales, the deuteron can be described within this pionless EFT. For example, one can 
calculate the deuteron form factors at momenta up to the pion mass [9]. This pionless 
EFT is limited in energy, but it is very simple (since all interactions among nucleons are 
of contact character) and high-order calculations can be carried out. Therefore, although 
its expansion parameter is not particularly small, high precision can be reached easily. 

In this very-low-energy regime the EFT of the two-nucleon problem is not much more 
than a reformulation of the analysis in Section 3. The full benefits of an EFT framework 
will, however, be evident when we consider the regime of momenta comparable to the 
pion mass in the next section. 
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4.1 Effective Lagrangian 

Nucleons with momenta much smaller than the pion mass are non-relativistic, and in this 
case, it is convenient to redefine the nucleon fields so as to eliminate the term proportional 
to mjy from the Lagrangian. In so doing, one obtains an infinite tower of operators 
proportional to powers of p/m^ « 1. This widely-used heavy-fermion formalism [75, 76], 
is nothing but a Galilean-covariant expression of the usual non-relativistic expansion. 
Since the non-relativistic EFT must match the relativistic theory for p ~ m^r, Lorentz 
invariance relates various terms in the tower of (p/mjv) fc - suppressed effective operators. 
Thus, one way to construct the effective Lagrangian is to write the most general rotational- 
invariant non-relativistic Lagrangian, then to relate parameters by imposing this matching 
condition, or "reparameterization" invariance [77]. Alternatively, we can simply write a 
relativistic Lagrangian and then take the non-relativistic limit. 

The most general Lagrangian involving two nucleon fields N, N and a photon field 
that is invariant under Lorentz, parity, time reversal and U(l) gauge symmetries is 

£n,pc = N^v-D + ^-^v-Df-D^ + lS^S^D^D"] (66) 

m N > 

where Kq(k,i) is the isoscalar (isovector) anomalous magnetic moment, and are 
the nucleon velocity and spin [v 11 = (1,0) and S* M = (0, a/2) in the nucleon rest frame], 
Dp = <9 M + ieQjyA^ is the electromagnetic covariant derivative, with Q N = (1 + r z )/2 the 
nucleon charge matrix, and F^ = d^A v — d v A lx . Here, as in the following Lagrangians, 
". . ." denote terms with more derivatives, which give rise to other nucleon properties such 
as polarizabilities. 

When we relax the restriction of parity invariance, we can write additional terms, such 

as 

C n,pv = — + ai r z )S,Nd u F^ + ..., (67) 



m 



N 



where (a x ) is the isoscalar (isovector) anapole moment of the nucleon. These terms 
were discussed in Refs. [33, 34]; they appear in PV electron scattering but not in the 
processes we focus on here. 

For the two-nucleon system, we need to consider contact terms with four nucleon fields. 
The simplest parity-conserving (PC) interactions are 

Cpc,nn = ~C S NNNN + 2C T NS ll NNS f ,N + ..., (68) 

where Cg, Ct are dimensional coupling constants first introduced in Ref. [3]. Their 
projections onto the two iViV S waves are 

Cos — C,s — 3Ct 

C ot = Cs + C T . (69) 

These parameters are related to the respective scattering lengths, while higher- derivative 
operators give rise to additional parameters, such as S-wave effective ranges and P-wave 
scattering volumes [8]. 
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For future reference, it is also useful to write down the first-quantized NN potential 
arising from £pc,nn- To order 0(Q), we have 

V^(q,p)=C s + C T ffi'ff2, (70) 

Similarly, we can construct PV two-nucleon contact interactions. A detailed derivation 
appears in Appendix A and leads at 0(Q) to 

C PVtNN = -^{-dN^N N^a -iD_N + C 1 N*iDLN NViV 

x 

-C x ie iik NhD^ajN N^a k N 

-C 2 N*NN*T 3 a -W-N + C 2 N j iD i _N NW^N 

-C 2 ie ijk NnD^N N^r 3 a k N 

-C 3 N^ a N N^r a a ■ iD-N + CsNW a iD i _N N^T a a i N 
-C 3 ie ijk NWSDXajN N^r a a k N 
-C 4 NW 3 N N^a ■ iD^N + C 4 NW 3 iDtN JVVJV 
-C 4 ie ijk N^T Z iD%UjN N^a k N 

-C 5 lab N^r a N N^r b a ■ iD-N + C h X ah N ] T a iD l _N JVWiV 
-C 5 labie ijk N^T a iD\ajN N^r b a k N 

-C 6 ie ab3 N*T a N N^r b a ■ iD + N) + . . . , (71) 

where we have introduced the short-hand notation 

NnD^N = (iD^N) 1 N±N^ (iD^N) . (72) 

(in momentum space, iD+ and iDt give rise to the difference and sum, respectively, of the 
initial and final nucleon momenta). The effects of the weak interaction are represented 
by the LECs Cj. We have normalized the operators to a scale A x = 4nF n ~ 1 GeV, as 
would appear natural in a pionful theory. One might then anticipate that the are of 
order GpA^, ~ 10~ 5 . In fact, as discussed in Section 5, naive dimensional analysis (NDA) 
suggests that these quantities have the magnitude C; ~ (A x /F 7T ) 2 g 7T , where g w = 3.8 x 1CT 8 
sets the scale for non-leptonic weak interactions. One may also attempt to predict such 
constants using models (see Section 6) and compare with the experimentally determined 
linear combinations discussed above. 

The 0(Q) Lagrangian Cpv,nn gives rise to the potential in Eq. (5), which generates 
energy-independent S-P wave mixing as discussed earlier. Higher-derivative PV operators 
lead both to energy-dependence in the S-P mixing amplitudes as well as mixing involving 
higher partial waves. Given the level of complexity already appearing at 0(Q), we will 
not consider these higher-order terms. 

4.2 Amplitudes 

In processes involving a single nucleon, amplitudes can be expanded in loops. The situ- 
ation is more subtle when two or more nucleons are present [3]. This is due to the fact 
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that intermediate states that differ from initial states only by nucleon kinetic energies 
receive Oim^/p) enhancements. A resummation then must be performed, leading, e.g., 
to nuclear bound states, and it is not immediately obvious that such resummations can 
be done while maintaining the derivative expansion necessary to retain predictive power 
order by order. The large values for the NN scattering lengths, however, provide justi- 
fication for such a procedure [8, 9]. Before considering PV effects, it is helpful to review 
what resummation technique yields for the case of the strong NN interaction. 

In lowest order, the S-wave NN interaction can be represented via a contact term 

T 0t = CM. (73) 

Including the rescattering corrections, the full T-matrix is found to be 

T t (k) = CM + CMGv{k)CM + . . . 



C 0i (n) 4tt 



1 - Coi(n)G (k) m N - mN c w(p) ~ V 



ik J 



(74) 



where 



G„(*) = to o G o(r >-) = / (gs ^A = -=£0. + ») (75) 

is the zero-range Green's function, which displays the large nucleon mass in the numerator. 
Identifying the scattering length via 

-- = Pt-,-» (76) 

and using the relation 

mt(k) = —Tiik) (77) 

47T 

connecting the scattering and T-matrices, we find 

mJk) = 1 1 , = (78) 

It is important to note here that since Oj is a physical quantity, it cannot depend on the 
scale parameter \x and this invariance is observed in Eq. (76), wherein the \x dependence of 
the Green's function is canceled by the corresponding scale dependence in —Ait / tunC^h) . 

We observe that the resummation is at this order completely equivalent to the unita- 
rization that lead to Eq. (26), and one can show similarly that in the NN system inclusion 
of higher-derivative operators reproduces higher powers of energy in the effective-range 
expansion [8, 9]. 

It is straightforward to generalize the above calculation to account for electromagnetic 
interactions. As shown in Ref. [78] (see also Ref. [79]) the unitarized pp scattering 
amplitude has the form 

(h\ mN Co^Cfa+ik)) e x P 2i<7 
mA ) ~ 4tt 1 - CMGc{k) ' (79j 
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where rj + (k) = Ma /2k, 



C 2 v (x) 



2nx 



2itx 



- 1 



(80) 



is the usual Sommerfeld factor, gq = argr(£ + 1 + ir] + (k)) is the Coulomb phase shift, and 
the free Green's function Go(k) has also been replaced by its Coulomb analog 



Gc{k) = J f27T)3 JSL-JL 



le 



(81) 



Remarkably, this integral can be performed analytically, yielding 



G c (k) 



m N 

47T 



H + m N a ( H(ir) + (k)) - log ■ 



c 



Here ( is defined in terms of the Euler constant via ( = 2rr — 7^ and 

H(x) = ip(x) + — log a;. 

The resultant scattering amplitude has the form 



(82) 



(83) 



m s (k) = 



C 2 ( V+ (k))e 2 ^ 



47T 



m N Cos(n) 



— ji — m^a 



H(irj + (k)) - log 



-c 



C 2 v (v + (k))e 



2iao 



1 

aos 



where we have defined, as before, 



Hv+(k)) - log 



47T 



mjv7ra 



C - ikc%( v+ (k)) 



(84) 



a 0s (n) MC 0s (n) 



and 



Hv+(k)) = ReH(ir] + (k)). 



(85) 



(86) 



The experimental scattering length ac s in the presence of the Coulomb interaction is 
defined via 



acs 



+ m N a ( log 



-c , 



in which case the scattering amplitude takes its traditional lowest-order form 

C 2 v ( V+ (k))e^ 



m s (k) 



-m N aH(ir} + (k)Y 



(87) 



(88) 



Of course, Eq. (88) requires that the Coulomb-corrected scattering length be different 
from its non-Coulomb partner, and comparison of the experimental pp scattering length 



—7.82 fm — with its nn analog 



—18.8 fm — is roughly consistent with 



Eq. (87) if a reasonable cutoff is chosen (e.g., jx ~ 1 GeV). 
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Having unitarized the strong scattering amplitude, we can now proceed analogously 
for its parity-violating analog. The lowest-order S-P mixing amplitude is 



T 0SP = W sp(fi). (89) 

Inclusion of S-wave rescattering effects while neglecting P-wave scattering and Coulomb 
contributions yields the result 

Tsp(k) = WospOj.) + WospMGoWCM + ...= (90) 



Writing Eq. (90) in the form 



m = ~T SP (U)-_ 

m N Coi{fJ.) 

= _ J _ Xt _ ik = \ t m t (k), (91) 

(Xi 

we identify the physical (/i-independent) S-P wave mixing amplitude via 

Similarly, including the Coulomb interaction, we find for the unitarized weak amplitude 

Tqsp = WospMCtiy+me*^ ^ \fpC?,( V+ (k))e^+^ 

(1 - C 0s {fi)Gc{k)) —^-m N aacsH(iri + (k)y 

where we have again neglected the P-wave scattering, and have identified 

Xsp ~ w (94) 

as the physical mixing parameter. 

Having obtained fully unitarized forms, we can now proceed to evaluate the helicity- 
correlated cross sections, finding, as before, at the very lowest energies, 

Al = stzfz = ^mfW) , 4Mf . (95) 

Somewhat more involved, of course, are processes involving more than two nucleons. 
Besides the inherent calculational difficulty, interesting new physics arises when three 
nucleons can overlap. When pions are integrated out of the theory, three- nucleon in- 
teractions become significant. In fact, it has been shown [11] that its strong running 
requires that the non-derivative contact three-body force be included at leading order 
in the EFT, together with the non- derivative contact two-body forces considered above. 
This three-nucleon force acts only on the S1/2 channel, and provides a mechanism for tri- 
ton saturation. The existence of one three-body parameter in leading order is the reason 
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behind the phenomenological Phillips line. Note that most three-nucleon channels are 
free of a three-nucleon force up to high order, and can therefore be predicted to high 
accuracy with two-nucleon input only [10]. Similar renormalization might also take place 
in the four-nucleon system. It remains to be seen whether the same phenomenon also 
enhances PV few-body forces. We defer a detailed treatment of A > 3 PV forces and 
related renormalization issues to a future study. 



For processes in which p ~ m^, it is no longer sufficient to integrate the pions out of the 
effective theory. Incorporation of the pion as an explicit degree of freedom requires use of 
consistent PV chiral Lagrangian, which we develop in this section. 

5.1 Effective Lagrangian 

Chiral perturbation theory (x?T) provides a systematic expansion of physical observables 
in powers of small momenta and pion mass for systems with at most one nucleon [58, 76]. 
The interactions obtained from %PT can be used to build four-nucleon operators arising 
from pion exchange, though care must be taken to avoid double-counting the effects of 
multi-pion exchange in both operators and wave functions (see below). In the approach 
we follow here, pionic effects are generally included non-perturbatively. Strong irN inter- 
actions are derivative in nature, and thus scale as powers of p/A x . As a result, one can 
include them while maintaining a systematic derivative expansion [57]. By contrast, weak 
nN interactions need not involve derivatives, but the small scale associated with hadronic 
weak interactions (gv) implies that one needs at most one weak vertex. In addition, ex- 
plicit chiral symmetry-breaking effects associated with the up- and down-quark masses 
also enter perturbatively, since m n « A x . To incorporate all these effects, we require 
the most general effective Lagrangian to a given order in p containing local interactions 
parameterized by a priori unknown low-energy constants (LECs). The corrections from 
quark masses and loops are then included order by order. 

We give here the basic ingredients to our discussion. (For a review, see Ref. [82].) 
The nucleon mass is much larger than the pion mass m^, so we continue to employ 
a heavy-nucleon field. The pion fields n a , a = 1,2, 3, enter through 



where F n = 92.4 MeV is the pion decay constant [83]. This quantity allows us to construct 
chiral vector and axial-vector currents given by 



5 EFT with explicit pions 




(96) 



A 




+ 0(tt 3 ), 



respectively. 
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Chirally-symmetric strong interaction pionic effects can be incorporated into the pio- 
nless Lagrangian by substituting — > V^, where the chiral covariant derivative is 

V li = D lt + (97) 

and by adding interactions involving A^. On the other hand, the quark mass matrix 
M. = diag(m u ,mrf) generates chiral-symmetry breaking that can be incorporated via 

X± = ?X?±&Z, (98) 

where 

X = 2B(s + ip), (99) 

with B a constant with dimensions of mass, and s,p representing scalar and pseudoscalar 
source fields. In the present application s = M. and p = 0, and in the following, we work 
in the isospin-symmetric limit, m u = = m. Isospin-breaking effects will generate small 
( < 1CT 2 ) multiplicative corrections to the tiny PV effects of interest here, so we safely 
neglect them. In this case, to leading order in the chiral expansion we have 

X+ = 2Bm + 0(n 2 ) 

X- = 2Bm— + 0^). (100) 

The building blocks for including a A field in the Lagrangian can be found in Ref. [84]. 
For simplicity we here integrate out A isobars. It is straightforward but tedious to use 
these building blocks to extend the results of our paper by including explicit A effects. 

We group terms in Lagrangians labeled by the chiral index v = d + f/2 — 2, where 
d is the number of derivatives and powers of the pion mass and / the number of fermion 
fields. We only display terms that are relevant for the arguments that follow. 

• Parity-conserving nN Lagrangian 
We then arrive at 

4°iv,pc = N[iv ■ V + 2g° A S ■ A]N , (101) 

with the lowest index. Similarly, we have for the next to leading order (NLO) 
Lagrangian 

C { Xpc = ^N{(v -V) 2 -V 2 + [S tl ,S u ][V^V} (102) 
- 2ig° A (S -VvA + vAS-V) + 2(k + mr 3 )e^ a pv a S^F^}N + ... 
and the next to next to leading order (NNLO) terms 

-.0 „o 



4\ ({S-V,vA} v-V + h.c.) — -^V (S-AV 2 + h.c.) 

2m 2 M Am 2 - v 



9°A 



2m 2 N 



S-V A-V + h.c.) + 2d w S ■ A( x+ ) 



+2d 17 (S ■ A X+ ) + td 18 [S ■ V, x-] + id 19 [S ■ V, ( x _)}} N + . . .(103) 
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Here the ellipses denote counter-terms not relevant in our present calculation, a 
complete list of which is given in Ref. [85]. The superscript "0" in g A and ji N 
indicates that these quantities must be appended by the corresponding loop contri- 
butions in order to obtain the physical (renormalized) axial coupling and nucleon 
magnetic moment. 

• Parity-conserving iViV Lagrangian 

For nuclear systems, we require the PC Lagrangian involving more than two nucleon 
fields. Here we will only need the lowest index [y = 0) terms, containing four nucleon 
fields. The relevant Lagrangian has the same form as Eq. (68), 

£% PC = -l-C s NNNN + 2C T NS^NNS^N + ..., (104) 

but here C$, Ct are constants whose numerical values are different from the ones 
in the pionless theory. This is because we are now removing soft-pion contributions 
from the counter-terms, and including them explicitly. 

The NLO four-nucleon corrections occur at v — 2, which will not be used since in this 
work we truncate the chiral expansion of the PV potential at 0(Q). Likewise, six- 
nucleon PC interactions first appear at v — 1 so their contribution to PV observables 
will be at higher order in loop diagrams. 

• Parity-violating ttN Lagrangian 

The lowest-index {y = —1) PV interaction arises from the irNN Yukawa interaction, 

_ KnN y3 at 

= -ih\ NN {pnix + - np-K~) + . . . (105) 

where 

= fVf - fr 3 ft (106) 

and the ". . ." denote the terms in this operator containing additional (odd) numbers 
of pions. At v — there exist also PV vector and axial-vector nNN interactions, 
detailed expressions for which can be found in Refs. [31, 32]. However, as discussed 
in Ref. [38], the effects of the vector operators can be eliminated through 0(Q) by 
using the equations of motion and by suitably re-defining the constant C n (defined 
below). The PV axial- vector couplings involve two or more pions, and, as pointed 
out in Ref. [32], such couplings renormalize h\ NN at 0(Q 3 ). Consequently, their 
contribution to the PV NN potential appears at 0(Q 2 ), via loop effects. 

At NNLO [y = 1) we find several new PV nNN operators that will contribute to 
the PV iViV potential at 0(Q). In principle, these operators can be expressed in 
terms of the quantities X£ R defined in Ref. [31], thereby allowing one to determine 
the full, non-linear dependence on the pion fields. For our purposes, however, it 
is sufficient to truncate the expansion of these operators at one power of the pion 
field, since terms containing additional pion fields only contribute to the PV NN 
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interaction beyond 0(Q). After implementing the strictures of reparameterization 
invariance, we obtain the Lagrangian 

QsN, pv = 2 ff L W^ M (r x *) 8 3 VS'JV + f^N [D"D X , (r x vr) 3 ] N 

X n X n 

+ ^hg^N(rxn) 3 N+... , (107) 
A x^n 

where we have chosen a normalization such that the constants k]^N ought to be of 
order a few times g n according to naive dimensional analysis (see below) and where 
the ". . ." indicate terms involving more than one pion field. 

Nominally, then, there exist three new, independent operators that contribute to 
the PV NN potential at 0(Q). A proof of their independence under reparameteri- 
zation invariance, following the arguments of Ref. [87], will appear in a forthcoming 
publication and we do not reproduce the full arguments here. Heuristically, how- 
ever, the existence of these operators can be seen from their correspondence with 
the independent 0(Q 2 ) scalars that can be formed from the independent momenta, 
nucleon spin, and pion mass 5 : p-ff, o - px p*, (p — ft) 2 , and m 2 . Naively, then, one 
would have expected four independent 0(Q 2 ) operators, rather than just three as 
given in Eq. (107), with the operator corresponding to p* • p given by 

ivS A (r x n) 3 l5 x N . (108) 

However, in a relativistic formulation of the theory, the corresponding operator can 
be rewritten in terms of N(f x n) 3 N and N D X D\, (r x 7r) 3 N through suitable 
integrations by parts and application of the LO equations of motion 6 ; consequently, 
it cannot exist as an independent operator in the heavy baryon formulation. Indeed, 

similar arguments eliminate an analogous term, N D S ■ AL) X N, from the parity 
conserving Lagrangian. In contrast, the remaining terms in C\ NN ^ PV cannot be 
eliminated in the relativistic theory via such arguments and, thus, must exist as 
independent terms in the non-relativistic case. 

We also note that in order for EFT with non-relativistic nucleon fields to match 
the fully relativistic theory, the coefficients k]^N m general receive contributions 
proportional to h\ NN that arise from a non-relativistic reduction of the LO PV 
nNN Yukawa interaction in addition to contributions that represent bona fide 0(Q 2 ) 
effects. This situation is analogous to what occurs for the 0{Q 2 ) nucleon magnetic 
moment operator, whose coefficient /xjv — Qn + &n receives a contribution (the 
Dirac term) that is dictated by relativity and that is proportional to the 0(Q) 
constant (Qn) and a genuine, a priori unknown 0(Q 2 ) contribution (the Pauli 
term) parameterized by the anomalous magnetic moment. In the present case, only 
kl a N b N receive contributions proportional to h\ NN as dictated by relativity: 

7 ia _ KvN N ^x^JI. i h lb — hl NN A x F n | . . 

k « NN ~ AV2 M 2 N ' k « NN ~ 8V2 M 2 N ' (1U9) 



5 The presence of the single pion field leads to a pscudoscalar interaction with the nucleon. 
6 We thank Vincenzo Cirigliano for discussions on this point. 
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where the ". . ." indicate the unconstrained 0(Q 2 ) contributions. 

In practical terms, only two of the operators in Eq. (107) are likely to be experi- 
mentally distinguishable. In momentum space, the second and third terms can be 
written as independent linear combinations of (p— p i, ) 2 + m^ and m 2 . The latter acts 
like a chiral correction to h\ NNl so to 0{Q) in the EFT, it cannot be resolved ex- 
perimentally. When inserted into the PV NN potential, the former cancels the pion 
propagator, leading effectively to an 0{Q) contact operator that is indistinguishable 
from the SR operator proportional to C 6 . In contrast, the effects of the remaining 
operator involving kl a NN cannot be absorbed into the LO 7r exchange potential or 
any of the short-range 0{Q) operators. Its contribution to the potential has been 
given in Eq. (12). 

• Parity-violating ^nN Lagrangian 

Finally, there exists also a contact n^yNN interaction at v — 1 [38], 

£%n,pv = -ie^rPa^F^ + H.c. (110) 

• Parity-violating NN Lagrangian and 'yirNN Lagrangian 

The v = 1 PV four-nucleon terms have the same form as in Eq. (71) but with 
the gauge covariant derivatives replaced by X^, the gauge and chiral covariant 
derivatives. In this case, the coefficients Cj, Ci will differ numerically from those ap- 
propriate to the pionless theory, since in the latter case, the effects of pion exchange 
are incorporated into the operator coefficients. 



5.2 Power counting 

Throughout this work we use power-counting arguments to guide us in the task of identi- 
fying the most significant contributions to PV observables. Power counting is carried out 
under an implicit assumption about the size of the couplings of the EFT. It is assumed 
that the couplings are neither particularly small nor particularly large compared with 
"naive dimensional analysis" (NDA) [86] , in which LECs scale with F n and A x as 

//2 

x (A X F W ) 2 x (gX , (HI) 
where d, p, f = 2k, k and n are positive integers and 

9w ~ 17? ~ 3,8 x ir§ ' (112) 

In the absence of weak interactions (n = 0), the LECs scale with a large mass scale 
as (A x )~ u , where v — d + f /2 — 2 is the chiral index defined earlier. Hence, one obtains 
the ordering of operators in Q/A x described earlier. At energies that are small compared 
with the mass of W and Z bosons, weak interactions have a strength given by the Fermi 
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constant Gp = 1.16639 x 10~ 5 GeV -2 . The effective operators they entail are proportional 
to (powers of) the Fermi constant times the square of a mass scale. A natural scale is 
the pion decay constant, so we assume that these operators have coefficients of order of 
G F F 2 ~ 1CT 7 . In Eq. (112), we use g n = 3.8 x 1CT 8 because this scale appears naturally 
in quark-model estimates as in Ref. [19]. 

Here we limit ourselves to n — 1. Up to two derivatives, then, we have one PV nNN 
Yukawa coupling h\ NNl three NNLO PV irNN couplings k]°^N, ten short- distance LECS 
Ci,Ci, and one additional independent PV LEC C n if we consider PV photo-reactions. 
As emphasized earlier, only five independent combinations of Cj and Ci are relevant to 
low-energy PV observables in few-body systems, while the effects of all but one of the 
NNLO PV nNN operators can be absorbed into other terms in the potential. In practice, 
then, the inclusion of pions leads to a total of eight independent LECs. From Eq. (Ill), 
the expected size of the relevant PV couplings is 

Knn ~ (^)&r (113) 
QA ~ (^) 2 g« (114) 



The most challenging part of the power counting is to order the strong-interaction 
effects. Here we count powers of Q, where as above Q denotes a small quantity such 
as the pion mass m n , an external momentum p, or the electric charge, e. For example, 
the strong nNN vertex is counted as 0(Q), the PV Yukawa vertex is O(Q ), the pion 
propagator is (D(Q~ 2 ), and the four-nucleon vertices proportional to Cs,t are also counted 
as O(Q ). 

In the one-nucleon system, a loop integral / d 4 k can be simply counted as 0(Q i ). If 
there are two or more nucleons, this naive counting breaks down. The reason is that within 
nuclei nucleons are nearly on-shell. Thus, instead of being 0(Q), the q° component of the 
pion four-momentum in the one-pion-exchange (OPE) diagram shown in Fig. 2 is 0(Q 2 ), 
since 

q»=p»-p^VjfA, (116) 



where i, f label initial and final states. This simply means that in first approximation 
OPE is static. Now consider the loop diagram generated by the exchange of two pions 
between two nucleons, and focus on the dq° integral, which is, schematically, 



/ 



dq° 



2tt E/2 + q°- c?/2m N + ie E/2 - q° - <f /2m N + ie \ (q ) 2 - q 2 - m 2 + ie J 
E-q*/m N + ie {q 2 + ml) + -"' (U7) 



where E ~ 0(p 2 /m N ) is the nucleon kinetic energy. The ". . ." are contributions from 
the pion poles, which scale according to naive power counting, and other small terms. 
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Figure 2: Parity- conserving iterated one-pion-exchange diagram. A solid (dashed) line 
represents a nucleon (pion). The dotted line indicates the cut line which picks out the 
two-nucleon intermediate state. 



Yet, the term shown explicitly, stemming from the nucleon pole, represents an 0{m^/Q) 
enhancement over naive counting. 

This enhancement is more general than the specific diagram considered above. It is 
present in any diagram that represents a time ordering displaying an intermediate state 
with nucleons only. Such an intermediate state differs from the initial state only by a 
difference of nucleon kinetic energies, which is small because of the heavy nature of the 
nucleons. This type of intermediate state already appeared in the pionless EFT, and led 
to the resummation (74), which is equivalent to unitarization of the potential, i.e. to the 
solution of the Schrodinger equation. 

To carry out the resummation in the presence of explicit pions, two approaches have 
been proposed, which differ in the treatment of pion effects relative to the contact interac- 
tions. In the simplest approach [89], pion interactions are assumed to be small compared 
to the non-derivative contact interactions, and only the latter are resummed. Unfortu- 
nately, this assumption does not converge for all NN channels at momenta of the order 
of the pion mass [90]. In the other approach [3], non- derivative contact interactions are 
assumed to be comparable to OPE, and both interactions are resummed. In its original 
form, Weinberg's approach was proposed as an expansion of the potential. This approach 
appears to be successful in accounting for a broad array of nuclear observables [2], but 
it, too, has problems: iteration of the chiral-symmetry-breaking piece of OPE leads to 
inconsistent renormalization [89, 42]. 

Progress has been made recently in the understanding of the power counting relevant 
for NN scattering at Q ~ m n [42]. If an expansion is made around the chiral limit, the 
aforementioned problems are in principle resolved, and one obtains an expansion that is 
both consistent and converges. More work is necessary to test the new power counting, 
but at this stage we can see the reason for its success. The iteration of OPE in the chiral 
limit, together with the non-derivative contact interactions, makes the NN amplitude 
numerically similar to Weinberg's original proposal. Therefore, while unnecessarily re- 
summing higher-order terms, Weinberg's power counting can still be used to organize the 
potential. 

With this scheme, we separate Feynman diagrams into two classes: two-particle re- 
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(a) 



(b) 



Figure 3: OPE diagram that contributes to the long-range part of the PV potential. The 
filled circle indicates the PV 71NN Yukawa coupling. 



ducible (2PR) and two-particle irreducible (2PI). Only 2PR diagrams lead to the anoma- 
lous enhancement factor after loop integration discussed above. The 2PI diagrams, in 
contrast, do not contain shallow poles, so they have the same power counting as the 
one-nucleon system. With this classification in hand, one can use effective field theory to 
organize the calculation order by order. The sum of 2PI diagrams yields the potential, 
which is the kernel for the Lippman-Schwinger (LS) equation. Through iterations the 
2PR diagrams are generated. Solving the LS equation, or equivalently the Schrodinger 
equation, one arrives at the amplitude from which scattering can be calculated, and whose 
poles are the NN bound states. 

In this work we will follow Weinberg's formalism and derive the PV NN potential up to 
0(Q). Only the 2PI PV diagrams are included in the PV potential. All 2PR diagrams can 
be generated when the PV potential is inserted in the LS equation. In practice, the PV 
potential is much smaller than the strong potential so it can be treated as a perturbation. 
One can treat it as a PV operator and calculate the PV matrix element using the wave 
function from LS equation with the strong potential. The connection with the expansion 
of Ref. [42] is easily made by further expanding in powers of m 2 . 



Using the above power counting we construct the PV potential, classifying terms according 
to their size. We truncate the chiral expansion of PV potential at 0(Q), although the 
procedure can be carried out to higher orders in similar fashion. The PC potential has 
been derived to 0(Q 4 ) in Ref. [4]. 

At 0(Q^ 1 ), the only contribution comes from OPE diagrams of Fig. 3, where the PV 
vertex is the LO Yukawa interaction and the strong vertex arises from the operator in Eq. 
(101). These diagrams give rise to a long-range potential, V^(k): 



5.3 The PV NN Potential 




■i[ti x r 2 J 3 



(gi + a 2 ) • k 

k 2 + nil 



(118) 



where the —1 subscript denotes the chiral index of the corresponding amplitude and where 
k = pi-p[= p' 2 - p 2 - 
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(a) (b) W 

Figure 4: Corrections to the long-range PV NN potential from insertions of (a,b) higher- 
order PC irNN terms, which are denoted by the unfilled circle, and (c) loops. 



Subleading corrections arise from several sources. First, there are corrections to the 
long-range potential from corrections at the PC vertex (see Fig. 4). As discussed in 
Appendix B, the corrections involving rfi6,i8,i9 amount to a renormalization of the bare 
coupling g A while the term containing d 17 does not contribute. The remaining terms in 
Eqs. (102,103) are proportional to Qa and do not introduce any new unknown constants 
into the PV potential. Since their contributions are discussed in Appendix B, we do not 
reproduce them here. 

Qualitatively new corrections arise at 0(Q) from long-, medium-, and short-range 
effects, Vf\ R: Vf\ R and V^ pi g R , respectively. The NNLO long-range contributions arise 
from inserting the operators in Eq. (107) in the OPE diagrams (see Fig. 5). As noted in 
above, the effects of the operators proportional to k l ^ N can be absorbed in the potential 
through a suitable redefinition of h\ NN and Cq. The momentum space form associated 
with the remaining operator is 



9Ak\ a NN 
A F 2 



Ti x r 2 



0"l 



fl X Pl&2 



Qi 



+ (!<-> 2) 



+ ..., (119) 



where Pi ($) is the initial (final) momentum of the ith nucleon, = — Pi, and the 
". . ." denote the 0(Q) contributions generated by corrections to the strong 7iNN vertex 
through NNLO [see Eqs. (102,103) ]. Taking the Fourier transform of Eq. (119) leads, 
after some algebra, to the co-ordinate space potential in Eq. (12). In a similar way, one 
may evaluate the contributions to V^\ R generated by order Q 3 contributions to the parity 
conserving irNN vertex. 

The short-range part arises from 

i) the PV NN contact interactions in Fig. 6 and 

ii) possible chiral corrections to PC NN operators Cs,t, as shown in Fig. 7. 

The contact interactions have exactly the same form as Eq. (5), so we do not reproduce 
the expression here. Of course, the values of the Cj, Q differ from those in the pionless 
theory, where they effectively account for the effects of low-energy pion exchanges. In 
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(a) (b) W 

Figure 5: Corrections to the long-range PV NN potential from insertions of (a,b) higher- 
order PV irNN terms, which are denoted by the circled filled circle, and (c) loops. 




35 




Figure 8: PV TPE triangle diagrams that contribute to the medium-range PV NN in- 
teraction at 0(Q). 




Figure 9: PV TPE crossed diagrams that contribute to the medium-range PV NN inter- 
action at 0(Q). 



principle, one would expect these couplings to be renormalized by ix loop effects, as in 
the case of h\ NN . As we show in Appendix C, however, such loop effects vanish to 
0(Q). Similarly, PV loop corrections to the leading-order PC operators — illustrated in 
Fig. 7 — generate no corrections to the short-range couplings at this order. 

The medium-range part arises from the two-pion-exchange (TPE) diagrams, in- 
cluding 

i) the triangle diagrams in Fig. 8, 

ii) the crossed diagrams in Fig. 9, and 
hi) the box diagrams in Fig. 10. 

The evaluation of these diagrams is somewhat involved, and we give a detailed dis- 
cussion in Appendix D. Here, however, we note a few salient features of the calculation. 
First, the explicit form of the TPE potential is linked to the definition of OPE and the 
procedure to subtract the iterated OPE from the box diagrams. The slanted box diagrams 
are meant here as a representation of the full box diagram with the iterated static OPE 
subtracted (according to the procedure explained in Appendix D). Relativistic corrections 
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Figure 10: PV TPE box diagrams that contribute to the medium-range PV NN interac- 
tion at 0(Q). 



(beyond those in OPE) appear at higher orders. Next, we regulate the loop integrals us- 
ing dimensional regularization. The regulator-dependence is removed by the appropriate 
counter-terms, which in general have the form given in Eq. (5). The remaining, finite 
parts of the integrals contain terms "regular" - or polynomial - in momenta and m % 
and "irregular", or nonanalytic, terms. The former are indistinguishable from operators 
appearing in Eq. (5) (and higher-order parts of the potential), whereas the latter are 
uniquely identified with the loop integrals. In principle, one may choose to retain ex- 
plicitly any portion of the regular terms and absorb the remainder into the short-range 
LECs appearing in Eq. (5). The meaning of the Ci,Ci is, thus, scheme-dependent. Here, 
we adopt a scheme in which all of the regular terms are absorbed into the corresponding 
Ci, Ci, leaving only the irregular contributions explicitly in V^ P j^ R : 




(120) 



where 




(121) 



and 
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(122) 



Thus, the PV TPE amplitudes produce contributions with the same spin-isospin struc- 
ture as the contact interactions C 2 ,Cq. In fact, since the regulator-dependent and reg- 
ular parts of the amplitudes can be absorbed into V^, we would not expect any new 
spin-isospin dependence to emerge from the divergent TPE amplitudes. The spatial- 
dependence of the finite, non-analytic part, however, is qualitatively different. We discuss 
this difference below. 

Finally, we observe that there is no PV three-nucleon force to 0(Q). In connecting a 
third nucleon via a pion-exchange interaction, one increases the order of a given diagram 
by the same amount as if one added an additional loop. Consequently, the ingredients 
given in Section 5.1 allow at 0(Q) only tree-level three-nucleon diagrams that involve 
the leading order PC vertices and the PV Yukawa coupling. However, these diagrams 
cancel against recoil terms in the iteration of the two-nucleon potential. In fact, the situ- 
ation here is analogous to the 0(Q 2 ) PC three-nucleon force, where a similar cancellation 
occurs[3, 5]. As a result, if one employs an energy-independent potential (as is usually 
more convenient in few-body calculations), one may omit these three-nucleon diagrams. 
Non-trivial three-nucleon PV effects should appear only at 0{Q 2 ), which is beyond the 
order of our truncation here. 

5.4 EFT PV Potential: Qualitative Features 

As shown above, the PV NN potential to 0(Q) is given by Eqs. (37,118,119,120). The 
corresponding coordinate-space V FY (r) can be obtained straightforwardly by taking the 
Fourier transform of these expressions. On the basis of the power counting, one would 
expect the OPE potential LR to dominate in those channels where it contributes, 
unless h\ NN is anomalously small compared with the NDA estimate in Eq. (113). This 
potential is, of course, not new [16]. Several contributions arise with chiral index v — 1. 
Although they are all formally of the same order in power-counting, their effects may 
nevertheless be distinct due to the different operator structures and spatial ranges. The 
SR potential has already been discussed extensively in the treatment of the pionless EFT. 
Qualitatively, the only impact on the SR potential of including the pion as an explicit 
degree is that the numerical values of relevant combinations of the Cj and Ci will differ 
for the theory with pions. 

The two-pion exchange contribution V^Xm also appears at 0(Q). The result in Eq. 
(120) appears to be the first analytic expression for the PV TPE potential that is model- 
independent and consistent with the symmetries of QCD. Although studies of PV TPE 
effects have appeared previously in the literature (see, e.g., Ref [18]), direct comparison 
with our treatment is difficult. First, we have not been able to find an analytic expression 
in the literature. Second, two terms in the PV TPE amplitudes that depend strongly 
on the cutoff would have appeared explicitly had we not used dimensional regularization. 
This regulator, or cutoff, dependence requires inclusion of short-range counter-terms in 
order to guarantee that physical observables are regulator-independent. In the analysis 
of Ref. [18], however, no mention is made of the counter-terms, and we suspect that 
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the corresponding TPE potential is not cutoff-independent. Third, the component of the 
TPE amplitude unique to the loop diagrams is determined by chiral symmetry, and it is 
notoriously difficult to maintain this symmetry without using (see Ref. [5] for an 
illustrative example in the parity conserving three-nucleon sector). The situation for PV 
interactions closely mirrors the developments in the PC TPE NN potential, whose first 
derivation in accordance with chiral symmetry was given within EFT [4] , and whose form 
was recently clearly identified in a phase-shift analysis of iVW data [91]. 
The PV TPE contributes two spin-isospin operators. One, 

6 = i[T 1 XT 2 ] 3 (a 1 + a 2 )-q, (123) 

appears also in LR , but C^ ,Loop {q) is not a simple Yukawa function. The structure 
is also the same as the h 1 ' term in the DDH potential, where it is usually neglected. The 
other spin-isospin structure, 

2 = l±li (a, x B 2 ) ■ q, (124) 

has the structure of the ft,* -term in the DDH potential. In Fig. 11 we plot the momentum- 
dependence of the coefficients of the operators 6 [Eq. (123)] and 2 [Eq. (124)] for V^\ R , 
in comparison with the corresponding components of V^ LR and the DDH potential using 
DDH best values from Table 1. 

As expected on the basis of power counting, the OPE potential gives the largest effect 
for q ~ m n . As q increases, the TPE potential grows and eventually overcomes OPE. This 
feature can be understood simply from the more singular nature of TPE: while LR 
scales as q' 1 at large q (or r~ 2 at small r), V±\ /[R scales as q 1 (or r~ 4 ). In comparison 
with isovector cu-exchange term in DDH, the 2 component of V^\ R has qualitatively 
similar behavior at low-g (up to an overall phase). The rise with q is more rapid, how- 
ever, indicating a longer effective range than for cu-exchange. As pointed out above, 
the Oq component at distances r^l/m^ is missing in DDH, while it is not particularly 
small in m R . This component will generate an additional energy-dependence in the 
same channels OPE contributes. Presumably, the conventional practice of neglecting the 
TPE component leads to inconsistency in the analysis of experiments that probe the 6 
operator at different scales. We see no theoretical justification to neglect TPE. 

It may, perhaps, be surprising that the TPE contributions to 02,6 become numerically 
non-negligible compared to the OPE effect at relatively low-momentum. For example, 
when q is of the order of typical Fermi momentum for nuclei (~ 200 MeV), the TPE 
contribution to Oq is roughly one third the OPE contribution. One may wonder, therefore, 
whether the EFT converges too slowly to justify truncation at 0(Q). One should keep in 
mind, however, that TPE effects always appear in tandem with short-range components 
of the same order and that the latter properly compensate for the most singular part of 
the TPE contribution. 

New long-range, single pion-exchange terms also arise at subleading order. The struc- 
ture of the operator associated with k\ a NN - shown in Eq. (12) - is distinct from those 
appearing in LR , Vl'\ AR and V/"^ as well as from the operators appearing in the 
DDH potential. Additional structures are induced by relativistic corrections to the PV 



39 



1 




_4 l i i i i i i i i i i i i i i i l 

100 200 300 400 500 600 700 800 

q (MeV) 



Figure 11: Components of the PV NN potential (in units of g n 10~ 4 MeV~ 2 ) as function 
of the momentum transferred (in MeV): OPE (thick solid line); Cq component of TPE 
(long-dash line); C<i component of TPE (short-dash line); Cq component of DDH (thin 
solid line); C 2 component of DDH (dotted line). 
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ttNN Yukawa interaction, which are neglected in the DDH approach (see Appendix B). 
A consistent power counting, however, requires that one include them along with the SR 
and MR operators. 

Finally, one might also worry that we have not included A isobar contributions explic- 
itly since kia —itln is comparable to m n . Indeed, in our treatment, A effects are implicit 
in the LECs. Had we kept the A as an explicit degree of freedom, it would contribute to 
the two-body PV NN interaction solely via loops. Because the PV 7riVA interaction ver- 
tices are of D-wave character, loops that contain this new PV interaction are generically 
two orders higher than the corresponding nN loops containing the PV Yukawa coupling. 
Similarly, there would also be A contributions to the renormalization of h\ appearing 
in hR . Since experiments are sensitive only to the renormalized Yukawa couplings, 
the treatment of A loops will only affect the interpretation of h\ and not its extraction 
from experiment (see the last article in Ref. [38]). The only new contributions from the 
A to the PV iViV interaction would be in the TPE potential where the A appears be- 
tween two PC 7i N A vertices 7 . There would also be three- nucleon diagrams that are the 
PV version of the leading PC three- nucleon force [5]. The calculation of these effects is 
straightforward, and they introduce no new , a priori unknown PV couplings. We leave 
the "improved" version of the PV EFT containing these effects for the future when it may 
be required by phenomenological considerations. 

5.5 Currents 

As discussed earlier, any experimental program aimed at determining the PV low-energy 
constants will likely include electromagnetic processes. In order to maintain gauge invari- 
ance, one must include the appropriate set of meson-exchange-current operators. Typi- 
cally in nuclear physics, one expresses the requirements of gauge invariance through the 
continuity equation 

V-J=[H,p\, (125) 

where = (p, J). For the long- and medium-range components of the potential, a mini- 
mal set of current operators satisfying Eq. (125) can be obtained by inserting the photon 
on all charged lines in one- and two-pion-exchange diagrams. The meson exchange current 
(MEC) operator corresponding to V^^, Fig. 12, is given in Ref. [21]. The operators 
associated with V^\ R and V P \ AR are more involved [see Fig. 13(a-d)]. In particular, con- 
struction of the MEC operator associated with V^\ AR is technically arduous, as one must 
evaluate a large number of Feynman diagrams — a task which goes beyond the scope 
of the present study. Thus, we defer a derivation of these MEC operators to a future 
publication. 

The foregoing set of MECs constitute a minimal, model-independent set required to 
ensure that Eq. (125) is satisfied. In addition, one may consider MECs that satisfy Eq. 
(125) independently from the terms in the potential. At 0(Q), we find that there exists 
one such MEC that is not determined from V PV by gauge invariance. This operator is 

7 In the two-nucleon PV interaction, these are diagrams analogous to those in Figs. 9, 10 but with the 
A substituted for a nucleon on the line without a filled circle. 
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Figure 12: Long-range PV meson-exchange currents in leading order. A wavy lines rep- 
resents a photon. 




(a) (b) (c) (d) (e) 

Figure 13: Corrections to PV meson-exchange currents: OPE from minimal substitution 
in the sub-leading (a) PC and (b) PV nNN vertices, (c) TPE, (d) short-range contribution 
from minimal substitution in the PV contact interaction, and (e) OPE from new 'jrrNN 
vertex. Not all ordering and topologies are displayed. 
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obtained by OPE with an insertion of the operator from Eq. (110) [Fig. 13(e)], leading 
to the momentum space two-body current 

V^QtvNnCtt 



J = 

mArA x F 7r 

and to Eq. (13) in co-ordinate space. 



,02-52 (<?1 + 52 ) X <?1 , n 
Ti =s h (1 <-> 



si 



(126) 



6 Short-Distance Archeology: Correspondence with 
DDH and Beyond 

In the ideal situation, a systematic EFT treatment would use experimental low energy 
measurements in order to determine the counter-terms X t , \®' 1,2 , pt, h\ NN) k\ a NN) and 
C n entirely from data. As emphasized earlier, there exists in principle a program of low- 
energy few-body measurements which will yield at least five linear combinations of these 
constants. Alternately, one would ultimately hope to gain a theoretical understanding 
of the values of these constants (and their linear combinations) probed by experiment. 
However, obtaining reliable theoretical predictions is complicated, since the PV AW inter- 
action involves a non-trivial interplay of weak and non-perturbative strong interactions. 
Indeed, carrying out a first-principles calculation of the PV LECs is not yet possible, 
since lattice QCD techniques are not yet sufficiently advanced to address this problem. 
Consequently, in order to say anything about the LECs beyond NDA estimates, theorists 
have of necessity relied on model approaches. In this section we illustrate how the PV 
LECs can in principle be estimated from details of the short-range dynamics. 

Before proceeding further, it is useful to comment on the correspondence with, and 
difference from, the conventional DDH formalism in the treatment of short-distance PV 
physics. 

(i) The EFT approach is systematic and model-independent. No assumption is made 
about the dynamics underlying the short-range interactions in the EFT, whereas the 
DDH formalism relies on a light pseudoscalar- and vector-meson-exchange picture 
as indicated in Fig. 1. 

(ii) The LECs Ci_5 have a straightforward correspondence with the DDH PV meson- 
nucleon couplings h^j 1 , h° p ' 1,2 . In the EFT framework, however, Ci_5 could deviate 
strongly from the DDH values, as we illustrate below. 

(iii) In terms of the DDH meson-exchange language we have the constraints 

(jDDH qDDH 

(JDDH = (JDDH = 1 + Xwi (127) 
qDDH (jDDH qDDH 

r<DDH = r<DDH = r*DDH = 1 + Xp 5 (128) 

o 3 o 4 <_/ 5 

where Xp,u> denotes the ratio between tensor and vector couplings of p, cu meson- 
nucleon interaction. In our EFT approach, however, C1-5 constitute five LECs 
whose values need not be related to Ci_5 as in the DDH picture. 
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(iv) The DDH parameter h'^ is generally discarded since its "best value" is tiny. In 
EFT, on the other hand, h'^ contributes to C 6 , but C 6 need not be small since it 
can receive a contribution, e.g., from ao meson exchange. Moreover, although the 
operator accompanying Cq has the same spin-isospin structure as PV pion exchange, 
these interactions have different ranges and may in principle be distinguished as long 
as a sufficient range of energies is probed. 



6.1 Resonance saturation 

One popular model approach — which we adopt here for purely illustrative purposes — 
assumes that the short-distance dynamics is governed by the exchange of light meson reso- 
nances. This "resonance saturation" approach has some theoretical justification from the 
standpoint of the large- N c expansion, where N c denotes the number of colors in QCD [92]. 
It is also supported by several phenomenological studies. It is well known, for example, 
that in the 0(Q 4 ) chiral Lagrangian describing pseudoscalar interactions, the low-energy 
constants are well-described by the exchange of heavy mesons [93]. In particular, the 
charge radius of the pion receives roughly a 7% long-distance loop contribution, while 
the remaining 93% is saturated by t-channel exchange of the p°. Similarly, in the baryon 
sector, dispersion-relation analyses of the isovector and isoscalar nucleon electromagnetic 
form factors indicate important contributions from the lightest vector mesons [94]. Fi- 
nally, the primary features of the NN PC potential seem to be well described in such 
a picture [95]. Thus, it seems reasonable to assume that low- lying-meson exchange may 
play an important role in the short-distance physics associated with the PV LECs. 

With these observations in mind, we invoke resonance saturation to arrive at illustrative 
estimates for the PV LECs. The relevant Feynman diagram is the same as in Fig. 1, 
where the exchanged bosons include all possible heavy mesons with appropriate quantum 
numbers. Here, parity violation enters through one of the meson-nucleon interaction 
vertices. While the DDH framework includes only the lowest-lying vector mesons to 
describe the short-distance PV NN interaction, we also consider the exchange of ao(980), 
cti(1260) and /i (1285), as well as the radial excitations of these systems. Of course, the PV 
LECs receive additional contributions from higher resonances, correlated meson exchange, 
etc. However, we limit our consideration to this set, as it already suffices to illustrate to 
what extent the short- distance PV NN interaction can differ from the predictions of the 
DDH model. 

In order to estimate specific values of LECs in the framework of the meson-exchange 
model we require the corresponding PC and PV meson-nucleon Lagrangians: 

• Vector-meson exchange 

The parity-conserving vector-meson-nucleon interaction Lagrangian reads 



£pNN — 9 p nnN 



£-uNN — 9ujNnN 



£$nn - 9<i>nnN 



i Xp ■ u 

. 2m N 

I X.U} . v 
I X<j> ■ V 

Ip + 7. «V9 

. 2m N 



r ■ p^N, (129) 
u}»N, (130) 
<^N. (131) 
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The parity-violating vector-meson-nucleon (VNN) interaction Lagrangian is given 
inRef. [19]: 



rPV 
*-pNN 



r PV 

r PV 

l ~<l>NN 



= Nl^ 5 



h 2 



h° p r ■ p, + h\p\ + -j= (3r 3 pJ -t- Pi ) 



2^6 

N(TXp„) 3 a^q ul5 N, 
N^u, [hi + hln] N, 
A^75</V K + h \A N - 



N 



h' 1 

2m n 



(132) 

(133) 
(134) 



Note that we have adopted the convention for 75 following Ref. [96], which is 
different from that used in Ref. [19]. 

To our knowledge, the following contributions to the PV NN short- distance inter- 
action have not been discussed elsewhere in the literature. 



a (980)-meson exchange 



£fcSw = 9a NNNr ■ a N, 



(135) 



£aoNN = h ao Ni-f 5 (r x a ) 3 N. 



ai(1260)-meson exchange 



(136) 



(137) 



Note that the structure Nia^ u q u (f • a^) j 5 N is analogous to the weak-electricity 
form factor of nuclear beta decay; it is parity conserving but CP violating and 
hence is not included. The PV Lagrangian is 



rPV 
t-aiNN 



+N 



<^< + « + ^(3r 3 a 1 [;-r.<) 



N 



2m jv 



hi 



hls-a' + h^ + ^iSnaSo r ■ < ) 



N, (138) 



where a\ is the neutral component of a\ meson. 
/i(1285)-meson exchange 

^/xAW = 9hNNNl»l5flN, 



(139) 



N. 



(140) 



In principle, one may also include in such a model exchange of the radial excitations 
of p, u, 0, ao, ai, f\ mesons is also allowed. 
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6.2 LECs with DDH framework and beyond 

With these couplings in hand, we can identify our predictions for the various low-energy 
constants. 

If we consider only vector-meson exchange a la DDH, we have 



a 



DDH 



QDDH 
QDDH 



= 1 + Xo, 2 = 1,2, 



1 + X P i = 3-5, 



QDDH 

C? DH = -Khl 

QDDH _a ,1 

C° DH = -AX 

C? DH = -AX 



w 

qDdh _ —K p g pNN h'p. 



qDDH = _^ h 2 

5 2^ p 



where we have defined 



However, within the context of resonance saturation, these LECs could also receive 
contributions from radial excitations of rho and omega mesons, and from do-, ai-, /i- 
meson exchange. Hence we have, more generally, 



(142) 



Cl,2 


r<DDH 
— °1,2 


1 ^iiadiai 
+ °1,2 
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fiDDH 
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C3-5 


riDDH 
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C3-5 


riDDH 
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| QRadial 
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C 6 


r<DDH 
- ^6 


_|_ QRadial 





where 



C/ 1 - -A fl g flNN h° fl , 
C2 1 = -h h g hNN h}- h , 
C3 1 = —A ai g ai NNh° 



- -Aa x g ai NN (h° ai + ^01) ' 

C4 1 = —IS. ai g aiNN h l ai1 

C? = -K^NN^+h^), 



r<ai _ A ai 



5 — ■^J^yaiNN'l'ai- 
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Ct 1 = -^g ai NN(h 2 ai +h 5 ai ), (143) 

Similar relations will hold for the radial excitations. 
6.3 Estimates for Ci-q, Ci 5 

As noted above, arriving at reliable theoretical predictions for the PV LECs, even within 
the context of a model framework, is a formidable task, and one which certainly goes 
beyond the scope of the present work. Nevertheless, it is useful to have in hand educated 
guesses for their magnitudes and signs, if for no other reason than to provide benchmarks 
for comparison with experiment. To that end, we quote below both expectations based 
on naive dimensional analysis and values obtained from correspondence with the DDH 
model. Future work could include, for example, computing the weak couplings entering 
Eq. (143), thereby providing model estimates for the departures of the Cj and C,- L from 
their NDA or DDH values. 

There exist various values for the parity-conserving couplings g P NN, Xpi 9u>nn, and x^ 
quoted in the literature [97, 98, 99]. Fortunately, the combination g P NN(l + Xp) takes 
roughly the same value in different approaches: g P NNi)- + Xp) ~ 21. Likewise various 
approaches consistently yield a very small value for Xu>- It is thus reasonable to use the 
values Xp — 6, Q p nn = 3 or Xp — 3.7, g P NN = 4.5. A word of caution is in order here. 
The proper accounting of chiral symmetry in multi-pion contributions might affect the 
extractions of strong couplings. For example, the effect of u exchange in AW scattering 
is significantly reduced when correct TPE is considered [91]. As we emphasized earlier, 
the estimates here should be considered to yield only an educated guess for the order of 
magnitude of the LECs. The theoretical uncertainty from this exchange model is much 
larger than the choice of Xp an d 9 p nn- Here, we simply use Xp — 3.7, g P NN = 4.5, 
Xu) = —0.12, and g^NN = 14 to make our best guess. Results are given in Table 2 and 
should be used with due caution. 

7 Conclusions 

In summary, we have performed a systematic study of the parity-nonconserving nucleon- 
nucleon potential, and have suggested ways by which the present confused experimental 
situation can be resolved. We have proposed breaking this program into two separate 
pieces: 

i) Since the low-energy parity-violating potential involves five S-P wave mixing am- 
plitudes, we have constructed a simple local effective potential in order to reliably 
extract such quantities from experiments involving only the AW, Nd, or Na sys- 
tems. We have also suggested the critical experiments that are needed in order to 
successfully complete this task and have given explicit formulas which will express 
the mixing amplitudes in terms of experimental observables. We have also suggested 
a two-phase experimental program, where phase one would include six (or possibly 
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LECs 



Naive Dimensional Analysis 



Best values 



Range 



Ci 
C 2 

c 2 
c 3 
c 3 
c 4 

c 5 

c 6 



±158 
±158 
±158 
±158 
±158 
±158 
±158 
±158 
±158 
±158 
±158 



32 
28 
17 
15 
63 
296 
95 
1 

-11 
-51 



-95 - 
-84 - 
13 - 
11 - 
-63 - 
-296 
-289 
- 
-13 - 
-61 - 



• 172 

• 151 
32 
28 

• 171 
+ 803 

• 520 
1 

► -8 
-28 



Table 2: Estimates of ranges and best values for PV coupling constants Ci_6, Ci_5 (in 
units of g n = 3.8 x 10~ 8 ). 



seven) measurements needed to test the consistency of the pionless EFT and phase 
two would involve additional measurements needed to determine the pion-related 
parameters if necessitated by the results of phase one. 

ii) A second important facet of this program is to confront the extracted phenomeno- 
logical potential with theoretical expectations. For this task, we have systematically 
constructed a parity- violating nucleon-nucleon potential V PV (f) within the frame- 
work of effective field theory using the Weinberg counting scheme up to the order 
0(Q). The correspondence with, and difference from, the conventional DDH poten- 
tial were discussed. 

In order for this scheme to come to fruition additional work is required on several fronts. 
Experimentally it is critical to complete the key experiments, resulting in a confirmed 
and reliable set of low-energy phenomenological parameters. Once these parameters are 
known, it is important to use them in order to analyze the heavier nuclear systems and 
resolve the various existing conflicts. Doing so could have important implications for the 
applicability of EFT to other electroweak processes in heavy nuclei, such as neutrinoless 
double /5-decay[81]. Future work is also needed to understand the relation between the 
underlying effective weak potential V PY (r) and the effective phenomenological parameters 
p t , A t , A* and should involve the best available nuclear wave functions. What should result 
from this program is the resolution of the presently confusing experimental situation and a 
reliable form for the parity- violating nuclear potential, which we hope will set the standard 
for future work in this field. 
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Appendix A: The PV NN Contact Lagrangian 

Since we employ the heavy-fermion formalism, one can build the most general PV op- 
erators by using heavy-baryon fields directly. This approach, however, yields redundant 
operators, which then have to be eliminated by imposing reparameterization invariance. 
Alternatively, we can obtain the relevant operators starting from the relativistic theory, 
then performing a non-relativistic expansion. We use ipN, 4>n for the relativistic nucleon 
field and N, to denote the nucleon field after non-relativistic reduction. In general, 
there exist twelve possible PV and CP conserving iViV-interaction terms up to one deriva- 
tive, which we write as 



Oi = -^-ip N lia flu q u il) N il) N lYl5*pN 




"X 



2 = -^TipNli^^q v ^N^Nr^^i)N 



O s = -^■ipNr a 'y fl ipNipNT a 'y fl 'y5^N 



° 3 = Tj' l l J Nr a ia l , u q u ilj N ilj N T a -f^ 5 ijj N 



ij2 — — 

° 2 = -T^^nIi^n^nt^^^n 




X 



(144) 



04 = -AilNTzia^lpNilN^lnlb^N 



04 = T^^Nrzl^N^N^l^l^N 
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o 5 

o 5 
o 6 

where X ab is defined in Eq 
that all the isospin violation arises from the weak interaction, thus neglecting isospin 
violation from up and down quark mass difference and electromagnetic interactions, since 
corrections from such effects are typically around a few percent and negligible for our 
purposes. The isospin content of the above terms is transparent: the 1 ■ l,r ■ r terms 
conserve isospin, the piece with X ab carries AI = 2, and all remaining pieces change 
isospin by one unit. 

In order to understand the constraints that relativity imposes, we consider a simple 
example — the expansion of 0\ and 0\. Up to 0(Q) we have 

Oi = fr^l-NhNNha ■ iD_N + NhiDiNNh^N 

-ie iik N^liD^(T j NN^l(T k N], (145) 



X 
X 

= j^e ab3 i! N T a ifj N i! N r b i'j 5 ilj N 

X 

= ^e ab3 ^ N T a ^^ N i! N r b i(T IIU q v 'y 5 'ilj N , 

. (6). In writing down these PV operators, we have assumed 



61 = -^rie ijk N jf liD i ,a j NN j 'la k N, 

A 6 

x 

Note that the two relativistic structures 0\ and 0\ together yield three distinct non- 
relativistic spacetime forms. However, only two linear combinations of these forms are 
independent according to the strictures of relativity. On the other hand, if we had started 
from the non-relativistic theory and tried to write the most general effective Lagrangian, 
we would have naively identified each of these three structures as being independent 
and would have mistakenly postulated three, rather than two, LECs. The requirements 
imposed on the independence of various non-relativistic operators which follow from con- 
sistency with the relativistic theory is known as reparameterization invariance. Physically, 
this invariance amounts to stating that the non-relativistic theory should not contain more 
physics (e.g., LECs) than the relativistic one. Analogous situations occur in heavy-quark 
EFT and in the non-relativistic expansion of the nucleon kinetic operator in heavy-nucleon 
EFT. 

Similar results follow for the operators O2-5 and O2-5 in that each set {Oi, Oi} gener- 
ates two independent combinations of non-relativistic operators and, thus, two indepen- 
dent LECs. On the other hand, after non-relativistic reduction, Oq and Oq yield exactly 
the same form up to 0(Q). Hence, these structures yield only one independent LEC 
in the non-relativistic theory even though there are two different LECs in the original, 
relativistic theory. The new LEC is a linear combination of g & and g 6 . 
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The full PV contact heavy- nucleoli Lagrangian at 0(Q) in the NN sector can then be 
written in the form given in Eq. (71), where the LECs Cj are related to the relativistic 
couplings Qi via 

C1-5 = , (146) 

2m n 

Cl-5 = &-5 + ^91-5 , (147) 

C 6 = £ 6 - ^S6- (148) 

We thus have a total of ten PV LECs describing PV short-distance NN physics. For 
the purpose of characterizing PV operators in the NN system through 0(Q), these ten 
constants are sufficient. 



Appendix B: Corrections to V™ LR 

There are subleading corrections to OPE that arise from Fig. 4, where the strong vertex 
comes from subleading PC operators in Eqs. (102-103). In fact, when the v — 1 operators 
from Eq. (102) are inserted in Fig. 4, the resulting PV potential is naively of O(Q ). 
However, with 

v-q = q = ?kpit „ (Q 2 /m N ), (149) 
2m at 

with i (/) denoting the initial (final) nucleon, we have 

t/ pv _ • 9aKnn f n xr 2 \ (gj - pg / )o t i ■ (pi f + pu) - (1 +-> 2) 

Vla ' LR " \V2m%F n \^2~), ' { } 

where q = p 2 f — T>2i = Pu — Pif- Thus, this contribution enters at 0(Q) and must be 
included for consistency. 

Similarly, the second and third operators from Eq. (103) are nominally v — 2 but lead 
to corrections that are 0(Q 3 ), since they contain two kinetic operators from v • D or v ■ A. 
However, the insertion of the first, fourth and fifth v — 2 operators from Eq. (103) in 
Fig. 4 lead to contributions at 0(Q). We list these terms below. The contribution from 
the first operator in Eq. (103) reads 



V N 

The contribution from the fourth operator in Eq. (103) is 



^ V LR = ^V-V LR . (151) 



t/ pv ■ 9aK nn { t 1 xt 2 \ (f u + pj f )a 1 -q + (1^2) 

V lc, LR * z^o „ o L x2 , _2 ' ^ LtiZ ) 



8V2m 2 N F n \ 2 y 3 -f + m^ 
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Finally, the contribution from the fifth operator in Eq. (103) reads 

rPv _ ■ 9aK nn fnXT 2 \ (a 1 -p lf )(q-p li ) + (a 1 -p li )(q-p lf ) + (1^2) 



Now consider the operators associated with di6-ig. In the isospin-symmetric limit, to 
leading order in the chiral expansion we have 

S-A{ X+ )~^S-A 

(S-A x+ )~0 
[S-V, X -]~£S-A 

1 IT 

[S-V,(x-)]~#S-A. (154) 

TV 

As a consequence, the LECs (ii6,is,i9 simply renormalize the bare 7rNN coupling constant 
g A at order 0(Q 2 ) while LEC dyj does not contribute. Up to the truncation order 0(Q) 
the corrections from these LECs to PV NN potential are automatically taken into account 
as long as we use the renormalized (or physical) g A in Eq. (118). 

Another possible source of corrections to the long-range PV potential is the insertion of 
subleading PV operators in Fig. 3. As pointed out in the Section 5, the PV vector operator 
does not contribute due to vector-current conservation. The axial- vector operator involves 
two pions and leads to loop corrections at 0(Q 2 ). The contribution from the remaining 
PV operator proportional to k^ a NN has been discussed extensively in the main body of the 
paper. 

Many chiral loops exist at this order, from self-energy and PC- and PV-vertex correc- 
tions. The chiral loops do yield contributions 0(Q). However, these effects are included 
in the renormalization of qa [82] and of h\ NN [32]. 

Appendix C: Loop Corrections to V^\ R 

The contact PV interactions in Fig. 6 appear at 0(Q). From a simple counting of the 
chiral order of vertices , propagators, and loops, it is clear that loop corrections to these 
PV LECs, shown in Fig. 14, first appear at 0(Q 3 ), which is beyond our truncation order. 

Potentially more important are the loop corrections to the contact PC interactions, 
where one vertex is the PV Yukawa coupling of the pion to the nucleon. The relevant 
Feynman diagrams are shown in Fig. 7. 

Take the Cs NN contact interaction as an example. For diagram (a-1), the amplitude 
is nominally 0(Q), and reads 

•M h l Cs ^9A f d D k i(S!-k) i % 

1 al ~ ^2 F n J (27r) D v 1 -(p' 1 + k)+tev 1 -( Pl + k)+iek 2 -ml + ie 

= -hlNNCs^S? J™ sds [du /^A 

x hi ; (155) 

[k 2 + svi ■ k + s(l — u)vi ■ p'i + usvi ■ pi + m^] 3 ' 
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Figure 14: Possible chiral corrections to PV iViV contact interactions. 



where s has the dimensions of mass, and where we have Wick-rotated to Euclidean mo- 
menta in the second line. From this form it is clear that %M a \ <x Si ■ v\ = 0. The same 
argument holds for diagrams (a-2)-(a-4). 
For diagram (b-1), the amplitude reads 

... hl C s V2g A f d D k i(S 2 -k) 

iM bl ~ h 



n 2 F T J (2tt) d v 2 ■ (p2 + k) +ievi ■ (pi + k) +iek 2 - ml + ie 
x ,,, . ... , , .„ * 5» =0, (156) 



[k 2 + svi ■ k + s(l — u)v ■ P2 + usv ■ p'i 



m 



where we have used the fact that v\ — v 2 — v — (1, 0) for low-energy iViV interaction. 
Similarly, (b-2)-(b-4) vanish at 0(Q). 

There remains a third class of diagrams, (c-l)-(c-4). These are 2PR diagrams and their 
amplitudes do not vanish at 0(Q). For example, the amplitude corresponding to diagram 
(c-1) reads 



... hl C s V2g A f d D k i(S 2 ■ k) 
iM cl ~ h n - — / ——^ — — 



2 F n J (2tt) d v 2 ■ (p' 2 -k)+iev 1 - (p[ + k) + ie k 2 - ml + ie ' 

However, only the 2PI part of these diagrams should be included. In other words, the 
contribution from the two-nucleon intermediate state should be subtracted from the ampli- 
tude. This can be done in old-fashioned time-ordered perturbation theory. Alternatively, 
we may use the following identity to accomplish the subtraction easily: 

\ = p + 2n5 (v-k). (157) 

—v-k + ie v ■ k + ie 
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The second term corresponds to the two-nucleon pole, while the first term is free of the 
infrared enhancement discussed earlier. After subtracting the two-nucleon-pole contribu- 
tion, the modified amplitude for diagram (c-1) becomes 

■M h i°s^9A f d D k i(S 2 -k) i i 

1 cl ~ 71 2 F n J (2ir) D v 2 -(k-p' 2 )+iev 1 -(p' 1 + k)+iek 2 -m 2 n + ie 

Similarly, the 2PI parts of diagrams (c-2)-(c-4) vanish. We see that diagrams (c-l)-(c-4) 
can be generated from the PC O(Q ) Cs,t contact potential and leading-order PV OPE 
potential by iteration in the LS equation. 

In summary, the chiral-loop corrections to the PV short-range potential occur at 0(Q 2 ) 
or higher. 



Appendix D: Derivation of V^ V MR 

Of course, a consistent calculation in EFT must include all loop diagrams present to a 
given order. In this appendix we give some details of the evaluation of the diagrams in 
Figs. 8, 9, 10. We use dimensional regularization for simplicity. 
Let us consider first the triangle diagrams in Fig. 8. 

• Flavor-conserving case 

In this case the initial and final state on each nucleon line are the same, i.e., a 
proton remains a proton and a neutron remains a neutron. Diagrams (c) and (d) 
are mirror diagrams of (a) and (b) in Fig. 8. We focus on (a) and (b). The sum of 
their amplitudes reads 

V2g A hj r d D k NjTfa ■ (2k - g) iViiV 2 LS 2 • (2k - q)N 2 

4F, 3 J (2ti) d v(p 2 + k)(k 2 - 
^9aK(^ ,i r , d D k 



' M(a)+(b) 4F3 J (2tx) d v(p 2 + k)(k 2 -ml)[(k-qf - m\ 



rl roo r 

16 / dx / dy 
Jo Jo J 



4F3 [Jo Jo *J (2n) D 
(vi-k)(S 2 -k) 
[k 2 — y 2 — m 2 — x(l — x)q 2 Y' 

+ 8(S 2 . q )l yiy J aM ,- ix) J — 
x ,,, .., .J IT^ }=°- (158) 



[k 2 — y 2 — m 2 — x(l — x)q 

After momentum integration the first term contains a factor v± • S 2 — 0. The second 
term vanishes due to the x integration, since the integrand is a total derivative. 

Flavor-changing case: n <-> p 
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In this case the sum of the amplitude for diagram (a) and (b) reads 

\/2g A h\ f d D k {pv\ ■ (2k — q)n) (nS 2 ■ qp) 



iM (a)+{b) 



4F3 J {2it) d v (p 2 + k){k 2 - ml)[{k - q) 2 - m\\ 

= -*^§^L(q)(pn)(nS 2 -qp), (159) 

x 7r 



where the function L(q) is defined as in Eq. (122). The sum of the amplitude for 
diagram (c) and (d) is, likewise, 

iM {c)+[d) = r^pLLiq) (np) (pS 2 ■ qn) . (160) 
Note that S 2 ■ q ~ —\oi • q- 

Summing the four diagrams and converting to momentum-space operators we get 

w^kr ^ ra ^ ' q ^ • (i6i) 

Clearly the sum of triangle diagrams has the same Lorentz, isospin structure as the 
Cq contact term. 

Consider now the crossed-box diagrams in Fig. 9. 

• Flavor-conserving case: pp — > pp, nn — > nn 

In Fig. 9, contributions from diagrams (c) and (d) are equal to (a) and (b). For 
pp — > pp the sum of (a) and (b) leads to 

* 4 ^2^)frl (p'Pi ■ ?, S?]p) (p ] Slp) . (162) 

Note that, for pp — > pp, initial particles are identical. The operator form will 
generate both (a), (b) and their mirror diagrams simultaneously. 

For the nn — > nn channel, there is an extra minus sign from PV Yukawa vertex. 

V2L(q)^-e ijk (ntgVn) (rJa k n) . (163) 

Combining both pp — > pp and nn nn channels, we get 

~ L{ r A2F eijk {N j r 3 q^N) (jVW) 

V-^ X 

- L( r g ^ h F eijk (W^'iV) (NW 3 a k N) . (164) 
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• Flavor changing case: n — > p,p — > n 
The sum of diagrams (a)-(d) leads to 

[ ~ 3L(q) + H{q)] wk tab " ( N ^ aN ) ( N ^ b ° ■ * N ) > ( 165 ) 

where H(q) is defined as in Eq. (122). 

Finally, we discuss the box diagrams. As in Appendix C, we have to subtract the 
contribution from the two-nucleon intermediate state. The corresponding time-ordered 
diagrams are shown in Fig. 10. After the subtraction the 2PR part, we find: 

• Flavor-conserving case: p — > p and n — > n 
For np — > np the sum of all diagrams leads to 

i4V2L(q)^i- (nt[5x • ?, Sf]n) . (166) 

x 

Note that, for pn — * pn, there is an extra minus sign from the PV Yukawa vertex, 

-iW2L(?)^ (pt[ 5l . g , 5 f]p) ( n t5j n ) . (167) 

x 

Combining both channels, we get 

» x 

- ^^ 6^ (jV+gW) (iVVa^iV) . (168) 
» x 7r 

• Flavor-changing case: n <-> p 

The sum of all diagrams leads to the same result as in the crossed-box case, 

^ [-3L(g) + //(g)] ||^ 3 (iWiV) (iVW ■ giV) . (169) 

In summary, the sum of one-loop, TPE diagrams is 
L (q) 9aK 



27=2 ( iVW -^) 



+ [-3L(?) + f|^ 3 (iWiV) (jvV* • giv) 

- v^L^fl^e^iVtgVw) (N^r 3 a k N) , (170) 
which leads to the medium-range potential (120). 
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Appendix E: Illustrative Estimates 



Having a form of the weak parity- violating potential V PV (r) it is, of course, essential to 
complete the process by connecting with the S-matrix — i.e., expressing the phenomeno- 
logical parameters Aj, p t defined in Eq. (36) in terms of the fundamental ones — Cj, Ci 
defined in Eq. (37). This is a major undertaking and should involve the latest and best 
NN wave functions such as Argonne V18. The work is underway, but it will be some 
time until this process is completed[101]. Even after this connection has been completed, 
the results will be numerical in form. However, it is very useful to have an analytic form 
by which to understand the basic physics of this transformation and by which to make 
simple numerical estimates. For this purpose we shall employ simple phenomenological 
NN wave functions, as described below. 

Examination of the scattering matrix Eq. (31) reveals that the parameters X s>t are 
associated with the short- distance component while p t contains contributions from the 
both (long-distance) pion exchange as well as short distance effects. In the former case, 
since the interaction is short ranged we can use this feature in order to simplify the 
analysis. Thus, we can determine the shift in the deuteron wavefunction associated with 
parity violation by demanding orthogonality with the 3 Si scattering state, which yields, 
using the simple asymptotic form of the bound state wavefunction[102],[103] 

Mr) =\l + Pt(° P + *„) • -iV + X t (a p - B n ) ■ -iV)l sff~e~ ir (171) 

V Ztt r 

where 7 2 /M = 2.23 MeV is the deuteron binding energy. Now the shift generated by 
V PV (r) is found to be[102],[103] 



8^ d (r) ~ J d 3 r'G{r,7)V PV {r")^ d {r') 




'y\r—f\ 

V PV (nMr') 



J dWV PV {7)Mr') (172) 




where the last step is permitted by the short range of V py (f). Comparing Eqs. (172) 
and (171) yields then the identification 



7 X tX t = J dh'iB, - a 2 ) ■ ?V PW {7)Ur')Xtio (173) 



where we have included the normalized isospin wave function £ since the potential involves 
Tx, r 2 . When operating on such an isosinglet np state the PV potential can be written as 

V PY (7) = A[(c 1 -3C 3 )(a 1 -a 2 )-HV/ m (r)+2/ m (r)--iV) 

x 

+ (Ci-3(%)(£ix£ 2 ).vy ro (r)l (174) 



57 



where f m (r) is the Yukawa form 



fm(r) = 



m 2 e mr 
4nr 



defined in Eq. (7). Using the identity 



(?i x a 2 )-(l + (Ti • (J 2) = ~ #2) 



(175) 



Eq. 173 becomes 



In 



X 



2M 4tt_ 2 r 00 , 3 

l^AfT^-^Wo ^ 

-2(3d - CJUr) 3 ^ + (3d - 3C 3 - d + d)^¥^ r) 



dr 



PT 4 1 2Mm 2 ( 2m(6C 3 -3C 3 -2C 1 + Ci) + 7(15C 3 - 3d - 5d + Ci) 
Xt 12 47rA| \ 



(7 + m) 2 



(176) 



or 



Mm 2 ( 2m{<oC 3 -3C 3 -2C l + d) + 7(15d - 3d - hC x + d) 
6ttA3 ^ 



(177) 



., , (7 + m) 2 , 

In order to determine the singlet parameter A" p , we must use the 1 S'o np-scattering 
wave function instead of the deuteron, but the procedure is similar, yielding[102],[103] 

d7(k)Xs = i-^&j dV(*i - a 2 ) ■ r^ PV (r>s (r')x s 6 (178) 

and we can proceed similarly. In this case the potential becomes 

2 



A 3 

x 



(d + C 3 + 4C 5 )(<? 1 - a 2 ) • (-iV/ m (r) + 2f m (r) • -iV) 



+ (d + d + 4d)(<?iX,r 2 )-V/ m (r) 
and Eq. (178) is found to have the form 



(179) 



48ttA3 3 

A, 



roc 

<?i - °2) 2 Xs / drr 3 { 
Jo 

#i 5o (r) 



x 2[C 1 + C 3 + 4C 5 ]/ m (r) 



+ [d + d + d + d + 4(d + d; 



1 ffim(r) 
1 dr 



= ~l2 Xt 



1 2Mm 2 
36 4ttA3 



1 



(k 2 + m 2 f 



x 



cos 5 s (4A; 2 (d + d + 4C 5 ) + (d + d + C 3 + d + 4(d + C 5 ))(k 2 + 3m 2 )] 



2m 



+ — sin 5 S [(d + d + 4d) K + 3F) + (d + d + d + d + 4(d + d)K 



(180) 
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which, in the limit as k — > 0, yields the predicted value for A™ p — 

K P = -lp faff^*) - 6 J!t A 3 {3[Ci + C 1 + Cs + C 3 + 4(C 5 + (%)] 

X 

- 2m< p [2Ci + Ci + 2^3 + ^3 + 4(2C 5 + C' 5 )]} (181) 
Similarly, we may identify 

Af = lim df(k) = 6?r ^ A3 {3[d + C 1 + C 2 + C 2 + C 3 + C 3 + C 4 + C 4 - 2(C 5 + C 5 )] 

- 2maf [2d + C 1 + 2C 2 + C 2 + 2C 3 + C 3 + 2C 4 + C 4 - 2(2C 5 + (%)] } 
A™ = ~— lim d T( k ) = - M A , |3[Ci + Ci - C 2 - C 2 + C 3 + C3 - C 4 - C4 - 2(C 5 + (%)] 



- 2maT[2C 1 + C x ~ 2C 2 -C 2 + 2C 3 + C 3 - 2C 4 -C 4 - 2(2C 5 + C 5 



In order to evaluate the spin-conserving amplitude p t , we shall assume dominance of the 
long range pion component. The shift in the deuteron wave function is given by 

but now with 8 

V&W = ^ff^i ~ + *0 ■ -^^(r) (185) 

Of course, the meson which is exchanged is the pion so the short range assumption which 
permitted the replacement in Eq. (172) is not valid and we must perform the integration 
exactly. This process is straightforward but tedious[100]. Nevertheless, we can get a rough 
estimate by making a "heavy pion" approximation, whereby we can identify the constant 
Pt via 

T • M 



(182) 



h PtXt ~ "'327 / rfV((?1 + ^ ' (OMOxrfo (186) 



which leads to [105] 



FT _ 1 47r /- . ^ n 2 Knn9*nn f°° , 3 dU(r) 1 
2n 327r 3 v2 JO dr m~ 

PT C 1 Kg wNN 7 + 2^71^ 
" V27 8Xt 967 71 (7 + m.) 2 (18?) 



8 Hcre we have used the identity 



(fi x f 2 ) = -i(fi - 7^)^(1 + n ■ f 2 ) (184) 



59 



We find then the prediction 



9*nn l + Zm^ ! 
Pt - 12V2n(l + m«) 2KNN ( } 

At this point it is useful to obtain rough numerical estimates. This can be done by 
use of the numerical estimates given in Table 2. To make things tractable, we shall use 
the best values given therein. Since we are after only rough estimates and since the 
best values assume the DDH relationship — Eq. (9) between the tilde- and non-tilde- 
quantities, we shall express our results in terms of only the non-tilde numbers — a future 
complete evaluation should include the full dependence. Of course, these predictions are 
only within a model, but they has the advantage of allowing connection with previous 
theoretical estimates. In this way, we obtain the predictions 

A t = [-0.092(^3 - 0.014d] m- 1 

Af = [-0.087(C 3 + 4C 5 ) - 0.037Ci] m' 1 

Af = [-0.087(C 3 + C 4 - 2C 5 ) - 0.037(d + C 2 )] m; 1 

Af = [-0.087(C 3 - C 4 - 2C 5 ) - 0.037(d - C 2 )\ m' 1 

Pt = OMQKm- 1 (189) 

so that, using the best values from Table 2 we estimate 



A t = 


-2.39 x 10" 




= -3.41 x 10~ 7 fm 


s 


-1.12 x 10" 


-7 -1 

m n -- 


= -1.60 x 10~ 7 fm 


Af = 


-3.58 x 10" 


-7 -1 


= -5.22 x 10' 7 fm 


Af = 


-2.97 x 10" 




= -4.33 x 10~ 7 fm 


Pt = 


1.50 x 10~ 7 


m~ l = 


2.14 x 10~ 7 fm 



(190) 
(191) 

Again we emphasize that in arriving at the foregoing expressions, we have used the DDH 
relationships between the Ci and C{. In the more general case, one should obtain expres- 
sions containing roughly the linear combinations given in Eqs. (46). A similar caveat 
applies to the expressions below. 

At this point we note, however, that Af is an order of magnitude larger than the 
experimentally determined number, Eq. (55). The problem here is not with the couplings 
but with an important piece of physics which has thus far been neglected — short distance 
effects. There are two issues here. One is that the deuteron and A^A^ wave functions 
should be modified at short distances from the simple asymptotic form used up until this 
point in order to account for finite size effects. The second is the well-known feature of the 
Jastrow correlations that suppress the nucleon-nucleon wave function at short distance. 

In order to deal approximately with the short distance properties of the deuteron wave 
function, we modify the exponential form to become constant inside the deuteron radius 
#[102], [103] 

^Tl { -e~ jR r < R 




±e" 7r r> R 

r 



(192) 
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where 

,7 exp7i? 



is the modified normalization factor and we use R=1.6 fm. For the NN wave function we 
use [102], [103] 

!. sin Jp 2 +p\r 
A — *= — r <r s 
Vp 2 +P 2 o r . (193) 
sin pr 1 e^ r „ . „ V 7 

pr -^-+ip r * 

where we choose r s = 2.73 fm and p Q r s = 1.5. The normalization constant A(p) is found 
by requiring continuity of the wave function and its first derivative at r = r s 



= \fp*+Po r * sin pr s -pa s cos pr s 
sin Jp 2 + p 2 r s pr s {l + ipa s ) 



As to the Jastrow correlations we multiply the wave function by the simple phenomeno- 
logical form [106] 

0(r) = 1 - ce~ dr \ with c = 0.6, d = 3 fm" 2 (195) 

With these modifications we find the much more reasonable values for the constants Af ,rip 
and X t 

Af = [-0.011(C , 3 + C , 4-2C 5 ) - 0.004(d + C 2 )] m; 1 
A™" = [-0.011(C , 3-C4 + 2C 5 )- 0.004(C 1 -C , 2 )]m; 1 
A7 = [-0.011(C 3 + 4C 5 ) - 0.004d] m; 1 
A t = [-0M9C 3 - 0. 0003d] m' 1 (196) 

Using the best values from Table 2 we find then the benchmark values 



Af = 


-4.2 x lO^m; 1 = 


-6.1 x 10~ f 


; fm 


\T = 


-3.6 x lO^m; 1 = 


-5.3 x 10~ f 


; fm 


ynp = 


-1.3 x lO^m; 1 = 


-1.9 x 10~ f 


; fm 


A t = 


-4.7 x lO^m- 1 = 


-6.7 x 10~ f 


; fm 



(197) 

Since p t is a long distance effect, we use the same value as calculated previously as our 
benchmark number 

p t = 1.50 x 10~ 7 m; 1 = 2.14 x 10~ 7 fm (198) 

Obviously the value of Af is now in much better agreement with the experimental 
value Eq. (55). Of course, our rough estimate is no substitute for a reliable state of the 
art wave function evaluation. This has been done recently by Carlson et al. and yields, 
using the Argonne V18 wavefunctions[104] 

Af = [-0.008(C 3 + C 4 - 2C 5 ) - 0.003(d + C 2 )\ m' 1 (199) 
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in reasonable agreement with the value calculated in Eq. (196). Similar efforts should 
be directed toward evaluation of the remaining parameters using the best modern wave 
functions. 

We end our brief discussion here, but clearly this was merely a simplistic model calcu- 
lation. It is important to complete this process by using the best contemporary nucleon- 
nucleon wave functions with the most general EFT potential developed above, in order 
to allow the best possible restrictions to be placed on the unknown counter-terms. 
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